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Summary 



A computer laboratory is presented for assisting 
students in understanding inventory systems, construct- 
ing mathematical inventory models, and finding optimal 
decision rules. Four computer programs and the principles 
guiding their preparation are discussed: 



ISL-1: Balancing Carrying, Shortage, and 
Replenishing Costs 



ISL-2: The Deterministic Reorder Point— Lot Size System 



ISL-3 : The Probabilistic Reorder Point— Lot Size System 
ISL-4 : A General Inventory Systems Simulation 



Four manuals describing the use of the programs are 
included. Each manual specifies how the unit costs, the 
demands for inventory, and a variety of policy decisions are 
presented to the computer. It gives illustrations of actual 
inputs and outputs, and also suggestions on model construction 
and/or optimization. Observations are made on the methodology 
of production of computer laboratory manuals. 



The use of the laboratory by students is discussed. 
Observations are made on adapting the programs to any computer 
system, on testing the laboratory, and on extension of the 
work to other fields. 
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1. Introduction 



]_ t ]_ inventory systems is a field in which carrying* 
SEortage, and replenishment costs can be controlled by 
m nirin g appropriate decisions (see references HJ» l J» 

[3l, [5], and [6]). This field has a wide range of 
theoretical and practical applications. It is particu- 
larly suitable for the study of model building and the 
derivation of optimal decision rules. Students of 
inventory systems can extend their learning experience 
to other fields in Operations Research, Management 
Science, Industrial Engineering, Economics, Business 
Administration, etc. 

1.2 Computer time sharing systems can now be used 
economically as an aid in the teaching process. They 
seem to be particularly useful in simulating an ideal 
laboratory environment for the study of inventory 
systems. This report describes computer programs and 
manuals for such an inventory systems laboratory. 

1.3 The laboratory has been designed to give the student 
tEree major options for each inventory system under study? 

A. Simulation of inventory fluctuation 

B. Long term averages 

C. Table of total costs 

He can use the options to check whether he understands 
the inventory system (Option A), whether the model he has 
constructed seems to be correct (Option B), and whether he 
hag found reasonably good decision rules (Option C). 
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2. Inventory Systems Computer Programs 



2.1 Inventory systems are characterized according to such 
factors as costs (carrying, shortage, and replenishing), 
demands (deterministic or probabilistic), leadtime (zero 
or non-zero), and policy (when to order and how much). 

By varying these and other factors, the number of different 
inventory systems that can be constructed is practically 
infinite. This report deals witn only four inventory 
systems. Each is described in detail in an appropriate 
manual in the appendices. 

2.2 The first system (iSL-l) is intended mainly to give the 
student an appreciation for the need to balance carrying, 
shortage, and replenishing costs. In the second system 

( ISL-2) demand is constant and leadtime is zero. This is 
an excellent system for learning how to build a model and 
how to determine optimal decision rules. 

The third system (lSL-3)» which is an extension of 
ISL-2, all ows demand to be both deterministic and probabilistic. 

It is a relatively complex system for which it is rather 
difficult to construct a model. The problem of finding 
optimal decision rules is even more difficult. 

The fourth system (ISL-4) is a rather general inventory 
system. It allows the study of systems with a variety of 
properties and policies. These include lost sales, probabilistic 
leadtime , prescribed variable demand, and the inventory policies 
( z ,q), (z,Z), and (t,Z). Systems with some of these properties 
are kn own to be extremely difficult to analyze . The simulation 
option of Program ISL-4 is therefore particularly useful for the 
study of such systems. 

2.3 A number of principles were used in preparing the programs 
for the four inventory systems. 

A. The output from the computer will consist only of requests 

for data and instructions, and the execution of the instructions. 

B. As many data and instructions as possible will be presented at 
one time. 

C. A special feature for checking correct results is not necessary. 
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Principle A essentially means that the program 
will not provide text material of any sort. Such 
| material should he in the manual. Principle B ensures 

j concentrated decisions by the student and economic 

utilization of the computer. Principle C has been used 
in the belief that the student is the better Judge of 
correct results. 

2.4 The mathematical formulation of the programs is 
based on reference [5l« Program ISL-1 is essentially 
Example 1-0, pp. 4-7. Program ISL-2 is related to 
Section 5-l> PP* 79-84. Program ISL-3 is mostly based 
on Section 13* 3 » PP* 246-252. Most of the mathematical 
formulation of Program ISL-4 is new. Portions of the 
program are based on Section 12-3, pp. 217-224. 

The main coding and implementation was done in the 
BA.SIC language [4] on the GS-265 Time Sharing System. 

At the end of each manual the full listing of the program 
is given. All illustrations in the manuals exemplify 
inputs and outputs on the GE-265 system. 

The programs were also coded and implemented in the 
( ’ FORTRAN IV language on the IBM-709^ Batch Processing 

System. The listings of the programs are given in the 
Appendix E. The outputs from these programs, using the 
data following each listing, are almost identical with 
the corresponding outputs in the manuals. 



3* Inventory Systems Laboratory Manuals 

3*1 We consider the preparation of manuals to be a most 
important task in using computers in the teaching process 
The manuals in Appendices A to D should be regarded only 
as a first attempt to exemplify the methodology which is 
involved. 
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| 3.2 Several principles have guided us in preparing the 

manuals • , _ . . 

' a. »rhe manual should only deal with the laboratory. 

I Instruction on inventory systems should be given in the 

I classroom and/or in a textbook. 

I b. The manual should be read before the laboratory. 

I The student should come to the laboratory prepared with 

I data for the computer. 

I C. The manual should contain numerous illustrations 

I which provide examples of the precise inputs to the 

| computer and actual copies of output from the computer. 

I 3.3 Bach of the manuals in Appendices A to D includes 

f tHe following sections : 

f A. Presentation of data 

B. Illustrations 4 . . 

[ c. Construction of model and/or optimization 

\ d. Extensions and problems for solution 

I E. Listing of the BASIC program 



Section A describes the parameters, the controllable 
I _ variables, and the available output options. It gives 

[( ) the names of all these variables in the program. Actual 

examples of data and/or inputs are also given. 

: Section B provides illustrations of actual computer 

| runs with data and/or inputs, and outputs. 



I In Section C a few suggestions are given on the 

I construction of the model and/or on the method of finding 

optimal decision rules. 

I Section D discusses some of the limitations of the 

I system and sug g ests possible extensions. All these 

| extensions require changes in the computer program. The 

j section also includes several problems for solution using 

| the program as it is. 

I 
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Section E gives the full listing of the couputer 
program as it is coded in BASIC and stored in the GE-265 
system. The listing is mainly intended for the teacher. 

3.4 The actual production of a manual is a time consuming 
jot which requires extensive planning. Much thought 
must be given to the illustrations, to ensure exemplifica- 
tion of all the features and shortcomings of the computer 
program. 

We used the editing capabilities of the Time Sharing 
System for the production of portions of the manuals. It 
was hoped thet the complete manuals could be stored in the 
computer so that changes could easily be made. However, our 
experience has shown that this method of production was 
rather time consuming and costly. 



4. Use of the Laboratory by Students 

4.1 F ull use of the laboratory will be made at The Johns 
Hopkins University in the academic year 1968-69 in the 
undergraduate course "Elementary Inventory Systems." This 
course was not offered in 1967-68. 

4.2 The laboratory was pretested by graduate and undergraduate 
students at Johns Hopkins. It was also used in the author's 
courses "Advanced Inventory Systems" in the Fall of 19^7. on 

a number of occasions. This preliminary use seem3 to indicate 
that the laboratory is an effective teaching device. It was 
gratifying to see the smile on a student's face when the 
long term averages which he predicted were actu all y printed 
by the computer. It was interesting to watch another student 
who claimed to have found the optimal decision rules of a 
system. He used the option in the program which gave the 
costs in the neighborhood of his solution but to his surprise, 
the solution was not optimal. Eventually he found the correct 
solution and he also smiled. 
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4.3 The following observation can be made regarding the 

maximum benefit that can be derived from the laboratory. 

A. Students should be assigned several problems for 
solution in advance. They should come to the laboratory 
with several sets of data and a full list of the results 
they expect from the computer. 

B. Students should work in groups of, say, 3 to 5 students. 
Each group should be supervised by a laboratory assistant. 

C. Whenever it is appropriate, students should submit a 
brief report documenting their work. 



5 . Conclusion 

5.1 The inventory system laboratory described in this 
report is ready for use by students on the GE-265 Time 
Sharing System and on the IBM-709 1 *- Batch Processing System. 

It can be adapted for use on any other computer system with 
relative ease. 

5.2 Several principles have been suggested to guide the 
preparation of the computer programs and the manuals for 
the laboratory. These principles, as well as the laboratory 
as a whole, should now be subjected to detailed testing. 

3.3 The methodology developed for the inventory systems 
laboratory can be readily applied to other fields, m 
particular, the following fields within Operations Research 
seem to be especially suitable : queueing theory, game theory, 
response surface analysis, stochastic allocation models, 
reliability theory, and replacement models. 
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1. Introduction 



The main purpose of the inventory system used in this laboratory 
(System l) is to introduce the student to the balancing of carrying costs, 
shortage costs, and replenishment costs, by trial and error methods In 
other laboratories (ISL-2, ISL-3. ISL-4) the inventory systems can be used • 
for construction of mathematical models and for evaluation of optimal 
decision rules. 
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2. Properties of System 1 

2.1 Demands 

The pl annin g horizon is composed of N periods. The demand D(l) at 
the beginning of period I is known. The student supplies the values of N, 
D(l), ... >D(N) . 

2.2 Unit Cost 6 

The cost of carrying inventory is Cl per unit per period. The cost of 
shortage is C2 per unit per period. The cost of replenishing is C3 per 
replenishment. The student supplies the values of Cl, C2, and C3. 

2.3 Decisions 

The program allows for two types of decisions: 

A. The initial Inventory which is available at the beginning of the firut 
period, B(l). There is no replenishing cost associated with this amount. 

B. When should replenishments be made and in what quantities. Up to N 
replenishments may be specified. A cost of C3 is associated with each 
replenishment. The decision variables are 1(1), R(l(l)), . . . ,I(M), R(i(m)) 
where I designates a period at the beginning of which there is a replenish 
ment of R(l) and where M is the number of replenishments. 
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2.4 Constraints 

It Is assumed that the demands during the planning horizon are exactly 
met by replenishments. That is, 

D(l) + ... + D(N) - R(I(1)) + ... + R(l(M)) 

Thus the initial Inventory B(l) is equal to the amount on hand at the end 
of the horizon. 

2.5 Objectives 

For the given demands and unit costs the student must determine by 
trial and error the decisions which will give a minimum total cost of the 
Inventory system. 

3 . Presentation of Data 

3.1 Periods: N 

The' planning horizon is composed of a specified number of periods. 
These periods may be days, weeks, months, etc. The number of periods is 
given in a data statement: 

910 DMA < N > 

For example , for a horizon of 12 months the data statement is : 

910 DMA 12 

3.2 Demands: D(l), ...,D(N) 

Jn each period I there may arise a demand of D(l). The demand is 
given in a data statement: 

920 DMA < D(l)» •• *»D(N) > 

For examp l e , the in tons for the 12 month horizon would appear as : 

920 DMA 10,20,20,30,20,30,0,0,40,30,20,20. 
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(Note that for each period there must correspond a number designating the j 

demand for that period, even if that demand happens to be zero.) j 

3.3 Costs: Cl, C2, C3 j 

Die car rying cost, Cl, is measured in dollars per unit per period. j 

For example: $.20 per ton per month. The shortage cost, C2, is measured 

i 

in dollars per unit per period. For example: $5«00 per ton per month. 

Die replenishing cost, C3» is measured in dollars. For example: $10.00 j 

for each replenishment. The inventory costs are placed in a data statement, j 

| 

thus: j 



930 DATA. < C1,C2,C3 > 

For example, 

930 DATA .20, 5«00, 10.00 



3.4 initial inventory: B(l) 

The program permits the inclusion of any initial inventory, B(l), 
left over from a previous planning horizon. This information is given in 



a data statement: 

940 DATA < B(l) > 

For example, for an initial inventory of 10 tons the data statement will be: 
9lfO DATA 10 



3.5 Replenishments; When - I, Hew much - R(l) 

As with the initial inventory, the program also permits control of 
replenishments. One has to furnish M periods, l(l) to l(M), and with each 
period, I, a corresponding replenishment, R(l). This appears in a data 

statement as: 

950 DATA < I(l)» R(I(1)) I(M), R(I(M)) > 
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For example , suppose N = 12 months, and replenishments are desired at 
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the beginning of the months: 2,4,9* H with the quantities: 40, 80, 

70, and 50 respectively. The data statement would then appear as: 

950 DATA 2,140,4,80,9*70*11*50. 

(When additional lines are needed to accommodate all the replenishment 
information, the lines from 951 through 959 may ®l so U3ed.) 

3.6 Summary 

Five data line s are needed. The first three constitute the parameters 
of the inventory problem: the length of the planning horizon, the demands, 
and the unit costs. The last two lines are the controllable variables: 
the amount on hand at the beginning of the planning horizon and the 
replenishments during the horizon. The general form of the data statement is: 

910 DATA. < N > 

920 DATA < D(1)»...*D(N) > 

930 DATA < C1,C2,C3 > 

940 DATA < B(l) > 

950 DATA < l(l)» R(I(1))»....I(M)» R(I(M)) > 

For example: 

910 DATA 12 , ^ ^ 

920 DATA 10,20,20,30,20,30,0,0,40,30,20,20 
930 DATA .20,5*00,10.00 
940 DATA 10 

950 DATA 2,40,4,80,9*70*11*50 

In Section 5.2 (p.12) the printouts resulting from these data are given. 

4. Examples 

Three examples are given in this section. In all the examples the 

pi arm-tug horizon is composed of 12 months. On January 1 the demand is 

10 tons, and at the beginning of the other 11 months the demands are: 
20,20,30,20,30,0,0,1(0,30,20, and 20. The total annual demand is thus 21(0 tons. 
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The un it costs are also the sane in all three examples. The 
carrying cost is $0.20 per ton per month, the shortage cost is $5.00 per 
ton per month, and the replenishing cost is $10.00 per replenishment. 



In al l the examples in this section it is assumed that replenishments 
occur at the same interval of time and that the quantity replenished is 
constant . 



4.1 Example 1 

Suppose the initial inventory on January 1 is zero and suppose that 
4 replenishments of 60 tons each are scheduled for January 1, April 1, 
July 1, and October 1. The data and results would be: 



\ 910 L)A f A 12 

(, 920 DATA lU* 20 * 20 * 30 * 20 * 30 * 0 * 0 * 40 * 30 * 20 * 20 

|(~") 930 DATA * 20 * 5 *U 0 * 10 . 00 

[ 940 DATA U 

j 950 DATA 1 * 60 * 4 * 60 * 7 * 60 * 1 0 * 60 

KUM 
vJ AIT* 
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Example 2 illustrates how costs can be changed by a suitable change 
in replenishments. In Example 1 the total cost is $15^.00 per year, whereas 
in Example 2 it is $150.00. The student should note that the increased 
yepleni shment costs were more than compensated for by the decreased carrying 
costs . 

It. 3 Example 3 



8 






I 







Suppose now that replenishments are the same as in example 2 but that 
the initial inventory on January 1 is 10 tons. The data and the results are now: 
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910 DATA 12 

920 DATA 10,20*20,30*20*30*0*0*40* 30*20*20 
930 DATA .20*5.00* 10.00 
940 DATA 10 

950 DATA 1,40*3* 40* 5*40* 7* 40* 9*40* 11* 40 
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0 


0 


40 


9 


40 


40 . 


40 


40 


10 


40 


0 


30 


10 


1 1 


10 


40 


20 


30 


12 


30 


0 


20 


10 




CARRYING 


SHORTAGE 


REPLEM I S8MEMT 


TOTAL 


AVERAGE 


25.8333 


0 


. 5 




LIMIT COST 


• 2 


5 


10 




COST /PEKl OD 


5.16667 


0 * 


5 


10. 1 


TOTAL COST 


62. 


0 


60 


122. 



TIMES 



1 SECS. 



/ 1 



awtttWKawaa***^ w 









vn 



( ) 



4.4 Summary 

The costs in the three examples can now he summarized. 



CARRYING 


SHORTAGE 


REPLENISHING 


TOTAL 


64 


50 


40 


i;a 


40 


50 


6o 


150 


62 


0 


6o 


122 







We note that by appropriate decisions any one of the three types of- costs 
can be increased (or decreased). This usually results in a decrease 
(or an increase) in some other cost. The problem of finding the best 
decision is discussed in the next section. 

In passing one should note that in every example the annual 
replenishments equalled ZkO tons - the total annual demand. If the 
totals were not equal, say. if one replenishment of 2 3 0 tons is scheduled, 
then the following result will he printed: 



910 DATA 12 

920 DATA 1 0* 20* 20# 30* 20* 30# 0* 0* 40* 30# 20# 20 
930 DATA *20# 5* 00* 10* 00 
940 DATA 10 
950 DATA 1*230 



KUM 
iiJA IT* 



ISL-1 



11:57 



03/04/67 



KEPLEM ISHMEM TS APE MOT EQJAL DEMAMDS 
CHECK LIMES 920 AMD 950 




1 



! 




TIME: 



1 SECS. 













! 



n 



i- 



! 



5 . Optimization 

Analytic methods of optimization are discussed in other manuals. The 
student should use program ISL-1 to find an optimal decision by trial and error. 
The examples in this section illustrate a trial and error approach. 



5.1 Example 4 

A careful study of the results of Example 3 indicates that inventories can 
be drastically reduced if replenishments are not constant but are equal to the 
corresponding 2 months* demands. Thus, for no initial inventory one can get: 



91tf DATA 12 

920 DATA 10# 20# 20# 30# 20# 30# 0# 0# 40# 30# 20# 20 
930 DATA . 20# 5. 00# 10. 00 
940 DATA 0 

950 DATA 1 #30# 3# 50# 5# 50# 9# 70# 1 1# 40 




KUM 
0 AIT* 



ISL-1 



12:02 



03/04/67 



PEKIOD 


BE G I M 


KEPLEM I Sri ME M T 


DEM AMD 


EM D 


1 


0 


30 


10 


20 


2 


20 


0 . 


20 


0 


3 


0 


50 


20 


30 


4 


30 


0 


30 


0 


5 


0 


50 


20 


30 


6 


30 


0 


30 


0 


7 


0 


0 


0 


0 


8 


0 


0 


0 


0 


9 


0 


70 


40 


30 


10 


30 


0 


30 


0 


1 1 


0 


40 


20 


20 


12 


20 

• 


0 


20 


0 




. CAKrtT IMG 


SHOKTAGE 


KErt-EM IS3MEM T 


• 

TOTi 


AVERAGE 


10*8333 


0 


. 416667 




LIMIT COST 


• 2 


5 


10 




COST /PEKI OD 


2* 16667 


0 


4* 16667 


6 • * 


TOTAL COST 


2 6 • 


0 


50* 


76. 



1 SECS. 


















1 



TTME: 



5.2 Example 5 



i 



Further trial and error attempts may eventually lead to the following 
results which seem to be optimal: 



910 DATA 12 

920 DATA 10* 20# 20# 30# 20# 30# 0# 0# 40* 30* 20* 20 
930 DATA .20* 5.00* 10.00 
940 DATA 10 

950 DATA 2*40*4*tSU*9*7U*ll*5U 



*UM 

dAIT. 



ISL-1 



12:09 



03/04/67 



PERIOD 


3EGIM 


rtEPLEM I SH MEM T 


DEM AMD 


EM D 


1 


10 


0 


10 


0 


2 


0 


40 


20 


20 


3 


20 


0 


20 


0 


4 


0 


30 


30 


50 


5 


50 


0 


20 


30 


6 


30 


0 


30 


0 


7 


0 


0 


0 


0 


3 


0 


0 


0 


0 


9 


0 


70 


40 


30 


10 


30 


0 


30 


0 


1 1 


0 


50 


20 


30 


12 


30 


0 


20 


10 




C Art* Y IMG 


SHORTAGE 


rtEPLEM I S4 ME Ml 


TOT, 


AVErtATE 


14. 1 667 


0 


. 333333 




JMIT COST 


.2 


5 


10 




COST/PEkIOD 


2.33333 


0 


3. 33333 


6. 


TOTAL COST 


34. 


0 


40. 


74 



TIME: 



1 SECS. 
















I 

y 





5.3 Example 6 

As another example consider an inventory system with a planning 
horizon of 7 days. Each weekday there is a demand for 5 box cars. On 
Saturday there is a demand for 10 cars. There is no demand on Sundays. 

The carrying cost is $10.00 per car per day. The shortage cost is $20.00 
per car per day. The replenishing cost is $1*00,00 per replenishment. 

What is an optimal replenishing policy? 

One can first attempt to balance carrying costs and replenishing 
costs while avoiding shortages. Suppose that only one replenishment is 
made per week. Then 35 cars would have to be ordered, say, on Monday. 

The carrying cost would then be $1000 per week and the replenisning cost 
would be $1*00 per week for a total of $1400 per week. 

For two replenishments per week, say 20 cars on Monday and 15 cars 
on Thursday, the carrying cost is $1*00 per week. The replenishing cost 
is $800 per week so that the total cost is $1200 per week. 

For thr ee replenishments per week the replenishing cost alone is 
$1200 per week. Hence, if shortages are to be avoided the optimal solution 
seems to be : 



910 DATA 7 

920 DATA 5, 5, 5, 5, 5, 1 0, U 
930 DATA 10,2 0,400 
940 DATA 0 
950 DATA 1,20, 5, 15 

KU.M 

WAIT. 



Jl I 






mm 



t*. ”, ^ 







08/04/ 67 



I5L-1 14803 



PE K I CD 

1 

2 

3 

4 

5 

6 
7 



BEGIM 

0 

15 

10 

5 

0 

10 

0 



CArtrfriMG 



AVEKAGE 5.71429 

UN I T COST 10 

COST /PERIOD 57.1429 



TOTAL COST 400. 



TIME? 1 SECS. 



kEPLENISHMEMT 

20 

0 

0 

0 

15 

0 

0 



SHORTAGE 

0 

20 

0 

0 



DEMAND 

5 

5 

5 

5 

5 

10 

0 



END 



15 

10 

5 

0 

10 

0 

0 



i 

1 



rtEPLEN I S 4 ME N 1 



.285714 

400 

1 1 4. 28 6 
800. 



TOTAL 



1 71. 429 

1 200. 



5.4 Example 7 

In Example 6 no shortages were allowed. One can now consider shortages 
as follows : 

Suppose that instead of two replenishments of 20 and 15 cars one chooses 
19 and 1 6 respectfully. This will reduce inventories on Mondays, Tuesdays, 
and Wednesdays, but will cause a shortage on Thursday. The net effect, 
however, will be a saving of 3 x 10 - 20 = $10 per week. Similar consider- 
ations lead to the following solution which seems to be the optimal solution. 



2 2 




nxa-iunw**.**' — 













910 DATA 
920 DATA 
930 DATA 
940 DATA 
950 DATA 


7 

5# 5# 5# 5# 5# 10*0 

10# 20# 400 
0 

1# 1 5# 5# 20 




206 

230 






KUM 

WAIT. 






ISL-1 


14:29 


03/04/67 




CAKKYIMG 


SHORTAGE 


AV/EKA3 E 


3. 57143 


. 714236 



dN I T COST 
COST /PEri ! OD 



10 

35.7 143 



20 

14. 2 >357 



KEPLEM ISHMEMT TOTAL 

• 2?357 1 4 
400 

114.236 164.236 



TOTAL COST 250. 



100 . 



300. 



1150. 



TIMES 1 SECS. 



Note that deletion of Lines 206 and 230 has eliminated the detailed 
printouts of the inventory fluctuations'. The student may wish to use this 
option when he is interested only in the summary of the results. 



6. Extensions and Problems for Solution 







6.1 Extensions 

System 1 is restricted in a number of ways. For example: 

A. Leadtime is assumed to be zero. 

B. Demands occur only at the beginning of each period. 

C. The cost of shortages are measured in dollars per unit per period. 

D. Demands are assumed to be known (deterministic) . 

It is possible to relax these restrictions. However, to do so one wo ul d 
need to change the Program ISL-1 (see Section 7). The student will find that 
changes in the properties A, B, and C above are relatively easy to handle. 
However, this is not the case with property D. 

6.2 Problems for Solution 

A. Find the solution of the system described in Section 4 (p. 5) when: 

Al. The carrying cost is $2.00 per ton per month (instead of $0.20 
per ton per month). 

A2. The shortage cost is $50.00 per ton per month (instead of $5.00 
per ton per month). 

A3» The replenishing cost is $100.00 per replenishment (instead of 
w $10.00 per replenishment). 

B. Find the solution of the system described in Example 6 (p. 13) when 
each unit cost is 10 times as large. (Solve these separate cases as Problem A 
above . ) 
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C. Demand is known to cycle over a l4-day period: it is 10 on days 1,2, 
and 3; 15 on days 4 to 7; 5 on days 8,9, and 10; and 7 on days 11 to l4. The 
carrying cost is $1.00 per unit per day, the shortage cost is $1.00 per unit 
per day, and the replenishing cost is $10. What is an optimal replenishing 

policy? 



D. Demand is uniform at the rate of 2400 points per year. The carrying 
cost is $0.56 per part per year. The ordering cost is $42 per order. No 
shortages are allowed. What is the optimal lot size? 



E. The weekly demands in an inventory system during one year are 



• 2,1,0,1,5,3.1»4,1,0,4,4,1,0,3»1»2,1,0,3»1»1»0»0,1»1»2,3»2,C,4,2,1,0,2,0,3*1»1* 



1,2,2, 3,0,1, 1,1,1, 3 • The carrying cost is $1.00 per pound per week, the shortage 
cost is $10.00 per pound per week, the replenishing cost is $30.00 per replenisn- 



ment. Find an optimal replenishing policy. 



F. Solve the system of Section 4 (p. 5) when the monthly demands are equal 
to 20 tons each. 



G. Solve Example 6 (p. 13) when the Saturday and Sunday demands are for 
5 cars (instead of 10 and 0 cars). 



X 4 ' 









7. Program ISL-1 



01 HEM I8L-1 F OH MAMJAL 1 

\d2 HEM BV E. nIADDO.a AMD O.BhAATEM# HE J0HM5 H0PXIM5 JM I VEH5I 1Y 
95 DIM BC 60 ># H( 60 ># 0( 60>#£C 6U> 



100 LE 1 0 = 0 
105 LET H = 0 
110 HEAD M 
115 FOH 1 = 1 TO M 
120 HEAD DC I ) 

122 LE f D = D + DC I ) 

125 MEX1 I 

127 HEAD Cl# C2# C3 

123 HEAD 8 C 1 ) 

130 HEAD I 

135 IF I =-99 THEM 150 
137 HEAD HC I ) 

140 LET H=h + hCI) 

142 30 TO 130 

150 IF D = H THE M 200 

155 PrtlMT "rtEPLEM 1 5rI'XE M 15 AkE MOT SUJAL DEMAMD8" 
160 PHIMT "CHECH LIME8 92u AMU 950" 



170 

200 
202 
204 
206 
208 
210 
220 
230 
240 
245 
247 
250 
270 
27 5 
230 
235 
290 
300 
305 
310 

315 

316 
320 
32 5 
330 
335 
340 
345 

350 

351 
910 
920 
930 
940 
950 
993 
999 



5TCP 

LET II = 0 
LET 12 = 0 
LET 13 = 0 

PitIMT "PEKIOD"# "BEGIM HEPLE M I SHMEM T"# "DEMAMD"# "EM D" 

PHIMT 

FOH I = 1 TO M 

LEI EC I ) = BCI) + HC I ) - DC I ) 

PHIMT I# BCI)# H( I )# DC I )#£C I ) 

IF EC I )> = 0 HEM 250 
LE T 12 = 12 - EC I) 

30 TO 270 

LET II = II + EC I ) 

IF H C I ) > 0 THEM 230 
30 TO 285 
LET 13 = 13 + 1 
LET BCI + 1 ) = EC I ) 

ME XT I 

LET II = Il/M 
LET 12 = 12/M 
LET 13=13 / M 
PHIMT 
PHIMT 

PHIMT " "# "CAkhTIMG SrIOKTAGE KEPLEM I SHMEM T"# "TOTAL 

PhIMT 

PHIMT "AVEKAGE"# II# 12# 13 
PhIMT "UMIT C05T"# C1#C2#C3 

PHIMT "C05T/PEHI0D"# I 1*C1# I2*C2# I3*C3# I HC 1 + I 2*C2+ I 3*C3 

PHIMT 

PHIMT "TOTAL C05T", IHCHM# I2*C2*M# I3^C3^M# 

PHIMT I HCl*tM+I2*C2*M+I3*C3*M 
DATA 12 

DATA 1 0, 20# 20# 30# 20# 30# 0# 0# 40# 30# 20# 20 
DATA .20# 5.00# 10. 00 
DATA 10 

DATA 2# 40# 4# 30# 9# 70# 1 1# 50 
DATA -99 
EM D 



2 << 
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1 . Introduction 



t This laboratory deals with a deterministic reorder point - 

| lot size system. The system is an extension of ISL-1 and forms 

j a basis for two other laboratories: ISL-3 and ISL-4. Program 

f ISL-2 has 3 distinct uses: 




1 



I 




A. For any reorder point and lot size the student can 
find out how the system behaves over a reasonable 
number of periods, and what the corresponding 
averages are . This is a way to determine whether 
he understands the system. 

B. The student can next attempt to build a model a. or 
the behavior of the system in the long run. lie 
can use the computer program to test his model by 
comparing numerical results that he derives on 
paper with the numerical results derived on the 
computer . 

C. Once the model of the system has been determined, 
the student can address himself to the problem of 
optimization: How to find the optimal values of 

the reorder point and the lot size. He must attempt 
•^0 find decision rules leading to the desired values. 
Any numerical results that he obtains can be checked 
with the third type of output available from the 
computer program. 




i 












! 










2. Presentation of Data 



The student, can control all the parameters in this program. 
Some will be fixed for several runs of the program, namely, the 
demand and the unit costs. These are entered in the primary input 



statement . 



The other parameters, which may vary with each run. 



are entered as secondary Input. They include such variables as 
the reorder point, the lot size, the initial inventory* etc. 



2.1 Demand and Unit Costs ('Primary Input ) 

The demand, U, can take on any positive value. Its dimension 
is in quantity units per time period, where the quantity units may 
be a ton, gallon, etc., and where the time period may be a day, 
a week, a month, etc. 



There are three types of unit costs. The unit carrying cost, 
Cl> is in dollars per unit per period. The unit shortage cost, C2, 
is also in dollars per unit per period. The unit replenishing cost, 
03, is in dollars per replenishment. 



The request for primary input is 
DEJ/AKD, CARRY, SHORT, RhPLfill 



in the form: 



? < U, Cl, 02, C3 > 



(The computer does not print the variable symbols, they are printed 
here as a visual aid to the student . ) 



Suppose the time period is one week and the quantity unit is a 
ton. Then, for example, we might have a demand of l j tons per week, a 
unit carrying cost of $1 per ton per week, a unit shortage cost of 



2*1 






jggsg 







v9 pox* ton per wool;, and a unit re plenishing cost of .,.36 per 
replenishment • Thus, the response to the* Input statement would bo; 

UJMAJID, CARRY , SHORT, Rh’PL 21 
? 5, 1, 9, 36 

2.2 Decision Variables and Options (Secondary Input) 

The program requires 10 secondary inputs during its running; 
the reorder point, S; the lot size, Q; the simulation index, XI; 
the initial inventory, Ql; the number of periods to be simulated, H, 
the number of periods to be printed, Nl; the long term averages 
index, X2; the total costs index, X3; the reorder point and lot 

size increment, Jl; and the primary input index, D. 

Uach of the indices XI, X2, X3> and D should be either 0 or 1. 

If XI » 1, then the simulation starts with an initial inventory 
of Ql units and is carried out over N periods. The details of the 
first Nl £ N periods are printed. In addition, a summary of averages 
for the N periods is also printed. If Nl = 0, only the summary is 
printed. 

When XI = 0, the simulation option is not used. In that case, 
the values of Ql, N, and III are irelevant; however, they must be 
supplied. 

If X2 = 1, then the long term averages summary is printed. 

This can be compared with the simulation summary. If X2 = 0, then 
the long term averages summary is not printed. 



4 








II' X3 - 1» then a table ol* costs is printed In the 
neighborhood, Jl» of the reorder point U und the lot uiao ti« 

If D - 1 , then at the end oi* the current run, the computer 

♦ 

will ask lor new primary input (:..e., new demand and unit coots 
will have to be supplied). Thereafter, or if D = 0, the computer 
will ask for new secondary input information. Thus, if the student 
wishes to retain the old primary input data and only change the 
secondary inputs on the next run of the program, he should enter 0 
as the value for D. 

The request for secondary input is in the form: 

t 

POIHT, LOT, SIM, HUT, HUH, PKEJVS, AVhSK, COSTS, ST JP, IATA 




\ 

! 






i 

! 



I 

j 

\ 

4 

1 



§ 





? < S, Q, XI, 01, H, HI, X£, X3, Jlf D > 

For example, the student may suppj'y the following inputs: 

POEIT, LOT, SIM, HOT, HUH, Pl^HJTS.AYJiR, COSTS, UTHP, DHTA 

7 "5, 20, 1, 10, 100, 8, 1, 1, *5, 0 

This input indicates that the student is interested in a reorder 
point of -5 and a lot sise of 20* He wants a simulation with an 
initial inventory of 10 units to be run over 100 periods. He asks 
that the details for the first 8 periods (and a summary for the 100 
periods) be printed. The student also asks for the long term averages 
summary. He also wants the table of total costs for reorder points 
ol -5*5, -5, and -4.5 units and for lot sises of 19* 9, 20, and 20.5 
units. Finally, he indicates that he would like another run with the 

same primary input. 
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3. Illustrations 
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The illustrations in this section are designed to exemplify 
all the capabilities of Program ISL-2 . They also form the basis 
for the next section which deals with models and their solutions. 

3.1 Balancing Carrying and Replenishing Costs 

Examples 3.11, 3*12, and 3*13 illustrate the inputs and 
outputs relating to an inventory system in which the demand is 
5 tons per month, the carrying ccst is $1 per ton per month, and 
the replenishing cost is $36. No shortages are allowed. Hence, 
in the primary input the shortage cost is chosen as $99999 

1 

per month. In all examples the reorder point is 5 tons. This 
ensures that no shortages occur. The lot size in all examples is 
20 tons. Therefore replenishments occur every k periods. 

In ibcample 3. 11 the secondary input calls for a simulation 
with an initial inventory of 15 tons. 10 periods are to be 
simulated and the details for all these periods are to be printed. 
The output is self-explanatory. Note that replenishments indeed 
occur whenever the reorder point is reached*. 

In iixample 3.12 there is a request for a simulation over 
100 periods and for the long term averages. The output is again 
self-explanatory. One observes that the simulation averages and 










and long term averages are identical. They dilfer, though, from 
the re cult s in Example 3»H* 

The secondary input of Sample 3*13 calls for total costs oi 
reorder points of 0, i>, and 10 tons, and lot sizes of 15, 20, and 
25 tons. The last entry in the input ensures a subsequent call 
for a new primary input, dote that the entry in the center of the 
table, :£4 per period, is the same as the cost obtained in 

sample 3*12. 



RUN 



ISL-2 



22S53 W2 SAT 12/23/67 

THE DETERMINI STIC REORDER POINT-LOT SIZE SYSTEM 



DEMAND# CARRY# SHORT# REPLEN 
? 5# 1# 99999# 36 



Illustration 3«1 



POINT# LOT SIZE# SIM# INI T# RUN# PRINTS# AVER# COSTS# STEP# DATA lacample 3*11 
j •>, 20# 1# lb# 10# 10' 0' 0' 0 ' 0 



REORDER POINT = 5 



LOT SIZE = 20 

SIMULATION 




PER* 


BEGIN 


DEMAND 


END 


CARRYING 


1 


15 


5 


10 


12* 5 


2 


10 


5 


5 


7* 5 


3 


25 


5 


20 


22* 5 


4 


20 


5 


15 


17* 5 


5 


15 


5 


10 


12* 5 


6 


10 


5 


5 


7. 5 


7 


25 


5 


20 


22* 5 


8 


20 


5 


15 


17* 5 


9 


15 


5 


10 


12* 5 


10 


10 


5 


5 


7. 5 



SHORTAGES 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



REPLENISHMENTS 

0 



1 

0 

0 

0 

1 

0 

0 

0 

1 



33 






£ 






















a simulation oh 



10 



PERIODS: 



7 



m.-ERAGES 
IN’T CO SIS 
COSTS PER PERIOD 



CARRYING 

1 A 
1 

1 A 



total costs per period = 2 a* a 



J temple 3.11 (cant’d) 



SHORTAGES REPLENISHMENTS 

0 . 3 

99999 36 

0 10.8 



+ ^ *>*"*> + *> 

r t f* v 4- *r t- 



^ 1* r 'f t t» r j* *h 7 f t 'f* »* 



POINT, LOT SIZE# 
? 5, 20, 



SIM, INI T# RUN, PRINT S, AVER# COSTS, STEP, DATA 

1 # 1 5, 1 00, 0, 1 , 0* 0^ 0 



Example 3.12 



REORDER POINT = 5 



LOT SIZE = 20 



SIMULATION 



FOR A SIMULATION OF 100 PERIODS: 



CARRYING SHORTAGES REPLENI SHMENTS 

AVERAGES 
UNIT COSTS 
COSTS PER PERIOD 

TOTAL COSTS PER PERIOD - 24 



THE LONG TERM AVERAGES 



0 . 25 

99999 36 

0 9 



CARRYING SHORTAGES REPL EN I SHMENTS 



AVERAGES 15 

UNIT COSTS 1 

COSTS PER PERIOD 15 



0 

99999 

0 



• 25 

36 

9 



TOTAL COSTS PER PERIOD = 24 

* * * * * * * * * * ** * * * * * * * $ # 



POINT, LOT SIZE, SIM, INI T, PIN, PRINTS, AVER, COSTS, STEP, DATA 
? b, 20, 0, 0, 0, 0, 0, i a c, j exa mple 3« 1 3 



REORDER POINT = 



LOT SIZE a 20 







5 






^«¥^!«rT>V(KO.W^S^ rr ^,^, r -_^~^-.--~ i ; i y _ 




H. POINT 
0 
5 

10 



0 T A L 


C 0 


ST T A b L 


E 




LOT 


SIZES 




15 




20 


25 


19. 5 




19 


19. 7 


2A. 5 




24 


24. 7 


29 • 5 




29 


29. 7 



iicsunple 3 . 13 (Coat'd) 



3.2 balancing Shortage and Replenishing Costs 

In Examples 3.21 and 3*22, as in the previous section, the 
demand is 5 tons per month, and the replenishing cost is £j6. 
However, this time no inventory is to be carried and the shortage 
cost is $9 per ton per month. The appropriate entries have been 
made in the primary input. 

The secondary inputs and the corresponding outputs in 
Examples 3.21 and 3.22 are self explanatory. Note the effect of 
the carrying cost on the total costs xn the table oi ..^cample 3» t -2< 



t r •»* 'i» r • 



DEMAND# CAKKY # SHOKT# KEPLEN 
? b# 99999# 9# 36 



Illustration 3*< 



POINT# LOT SI2E# SIM# I Nil# KUN# PKINTS# AVEK# COSTS# STEP# DATA 
? -10# 10# 1* 0* 10* A# 0# 0* 0* 0 



Example 3«21 



KEOKDEK POINT = -10 



LOT SIZE = 10 



SIMULATION 



PEK. 

1 

2 

3 

4 



BEGIN 
0 
- 5 
0 

-5 



DEMAND 

5 

5 

5 

b 



END 
-5 
-10 
-b 
- 10 



CAKKY IN G 
0 
0 
0 
0 



SHOKTAGfc.S 
2* b 
7. b 
2* b 
7. b 



KEPLEN I SHMENTS 
0 



1 

0 

1 




rr 
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FOR A SIMULATION OF 10 



PERIODS! 



CARRYING 



acaniple 3,21 (Cant'd) 

SHORTAGES REPLENI SHMENTS 



AVERAGES 
UNIT COSTS 
COSTS PER PERIOD 



0 

99999 

0 



5 

9 

AS 



• 5 
36 
18 



TOTAL COSTS PER PERIOD a 63 



POINT# LOT SIZE# SIM# INI T# RUN# PRINTS# AVER# COSTS# STEP# DATA 
? -10# 10# 1# 0# 100# 0# 1# l# 2*5# 1 Example 3»22 



REORDER POINT * -10 



LOT SIZE » 10 

SIMULATION 



FOR A SIMULATION OF 100 PERIODS! 



CARRYING 



SHORTAGES REPLENI SHMENTS 



AVERAGES 




0 


5 


• 5 


UNIT COSTS 




99999 


9 


36 


COSTS PER PERIOD 




0 


45 


18 


TOTAL COSTS PER PERIOD • 


63 








THE L 0 


N G 


TERM 


AVERAGES 






CARRYING 


SHORTAGES REPLENI SHMENTS 


AVERAGES 




0 


5 


• 5 


UNIT COSTS 




99999 


9 


36 


COSTS PER PERIOD 




0 


45 


18 


TOTAL COSTS PER PERIOD * 


63 








THE TOTAL 


C 0 


ST T A 


B L E 






LOT 


SIZES 






7. 5 




10 


12.5 




ft. POINT 

-12.5 170.25 




175. 5 


70.65 




-10 147.75 




63 


25050.2 




-7.5 57.75 




31293 


100034. 





‘A 



1 





TTTZ 
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3.3 balancing carrying and Shortage Costs 



,£camples 3. 31 and 3.32 illustrate inventory systems in which 



only the reorder point is subject to control. This is the reason 
for selecting a zero cost of replenishing in the primary input, 
and for equating the demand and the lot size (5 tons). 



The reader should attempt to explain the carrying and 
shortage averages. Otherwise, tue outputs are sell —explanatory. 



¥ *r i* *•• ♦ / ¥*fT t 

DEMAND# CARRY# SHORT# REPLEN 
? 5# 1# 9# 0 



+ * + + + m + + + + m 



Illustration 3*3 



POINT#LOT SI2E# SIM#INI I# KUN, PRINT S# AVER# CO STS# STEP# DATA 
? -2# 5# 1# 3# 10# 4# 0# 0# 0# 0 



l&ample 3.31 



REORDER POINT = -2 



LOT SIZE = 5 



SIMULATION 




BEGIN 

3 

3 

3 

3 



DEMAND 

5 

5 

5 

5 



END 

-2 

-2 

-2 

-2 




CARRYING 

.9 
.9 
• 9 
.9 



SHORTAGES 



• 4 

• A 

• A 

• A 



REHL ENI SHMENTS 
1 
1 
1 
1 



FOR A SIMULATION OF 10 



PERIODS: 



CARRYING 



SHORTAGES 



REPLENISHMENTS 



AVERAGES 
UNIT COSTS 
COSTS PER PERIOD 



.9 

1 

.9 



• A 
9 

3. 6 



1 

0 

0 



TOTAL COSTS PER PERIOD = 4.5 




**4#* '' ****** + 4r *> 

-57 
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POINT# LO T SIZE* SIM# INI 1# RUN# PRINTS# AVER# COSTS# STEP# DATA Example 3.32 
? “ 1# S# 1# 4# 1 0* 0# 1# 1# 0* 5# 1 



REORDER POINT = -1 



LOT SIZE = 5 

SIMULATION 



FOU A SIMULATION OF 10 



PERIODS: 



CARRYING 



SHORTAGES REH. EN I SHM EN T S 



AVERAGES 
UNIT COSTS 
COSTS PER PERIOD 



1*6 

1 

1.6 



• 1 

9 

.9 



1 

0 

0 



TOTAL COSTS PER PERIOD * 2.5 



THE LONG 



TERM 



AVERAGES 



CARRYING 



SHORTAGES REPLENISHMENTS 



AVERAGES 
UNIT COSTS 
COSTS PER PERIOD 



1.6 

1 

1.6 



• 1 

9 

• 9 



1 

0 

0 



TOTAL COSTS PER PERIOD = 2.5 



THE 


TOTAL 


C 0 


ST TABLE 






LOT 


SIZES 




4 . 5 




5 5. 5 


.POINT 






• 


1.5 


3. 25 




3.25 3.29 54(5 


1 


2. 361 1 1 




2.5 2.659 09 


. 5 


2. 02778 




2.25 2.47727 



^ j0* ^ 4# 4# <4r <*• ^ <1* - V *#* a* 4- * ■ * 
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3.4 Balanciiig Carrying, Shortage, and Replenishing Costs 

In a certain sens#, llcamples 3.41 and 3.42 are extensions 



of all previous examples. ho special constraints are 
imposed on the inventory system. The demand is 5 tons per 
month, the carrying cost is $1 per ton per month, the shortage 
cost is $9 per ton per month, and the replenishing cost is $3 

Previous examples seem to indicate that the optimal 
reorder point is about -1 ton and that the optimal lot size is 
20 tons. However, the table of juxample 3*42 shows that 
that optimum is elsewhere. 



DEMAND# CARRY # SHO RT# REPL EN 
? 5# 1# 9# 36 



Illustration 3*4 



POINT# LOT SIdE# SIM#INIT#RUN# PRINT S# AVEK# CO ST S# STEP# DATA 



? - I 



20 # 



1 # 



14 # 10 # 



10 # 



0 # 



0 # 



0 # 



0 



ibcample 3.4l 



REORDER POINT = - 1 



LOT SI2E « 20 

SIMULATION 




PEK. 


BEGIN 


DEMAND 


END 


CAK HYING 


SHORTAGES 


REPLENISHMENTS 


1 


14 


5 


9 


1 1. 5 


0 


0 


2 


9 


5 


4 


6. 5 


0 


0 


3 


A 


5 


- 1 


1*6 


• 1 


1 


4 


19 


5 


14 


16. 5 


0 


0 


5 


14 


5 


9 


1 1* 5 


0 


0 


6 


9 


5 


4 


6. 5 


0 


0 


7 


A 


5 


- 1 


1.6 


• 1 


1 


8 


19 


5 


14 


1 6* S 


0 


0 


9 


1 4 


5 


9 


11. S 


0 


0 


10 


.9 


5 


4 


6. 5 


0 


0 



<*? 



















FOR A SIMULATION OF 10 



PERIODS* 



Example 3.41 (Cont'd) 






AVERAGES 
| UNIT COSTS 
| COSTS PER PERIOD 



CARRYING 


SHORTAGES 


REPLENISHMENTS 


9.02 


• 02 


• 2 


1 


9 


36 


9.02 


• 18 


7.2 



TOTAL COSTS PER PERIOD = 1 6.4 



j 

I ************ ««*****$** 

I 

[ 

I 

POINT# LOT SIZE# SIM# INI T# RUN# PRINTS# AVER# COSTS# STEP# DATA Example 3. 1 # 
( ? - 1 # 20 # 0 # 0 # 0 # 0 # 1 # 1 ' l» 1 




REORDER POINT = -1 LOT SIZE * 20 

THE LONG TERM AVERAGES 





CARRYING 


SHORTAGES 


REPLENISHMENTS 


AVERAGES 


9.025 


• 025 


• 25 


UNIT COSTS 


1 


9 


36 


COSTS PER PERIOD 


9.025 


• 225 


9 



TOTAL COSTS PER PERIOD = 18.25 



THE 


TOTAL 


C 0 


ST T A B L 


E 






LOT 


SIZES 






19 




20 


21 


R. POINT 
-2 


18.0263 




18. 


18.0238 


-1 


18.2368 




18*25 


18. 309 5 


0 


18.9737 




19 


19.0714 







* 



3*5 Special Features 

In all illustrations considered thus far the inputs were 
such tiiat the reorder point was always reached precisely. That is, 
we never had a case where the amount on hand at the end of a period 
was Delow the reorder point, Example 3«i>l now illustrates such a 
case. Note also that the simulation results do not agree with the 
long term averages. 

Example 3*51 also shows that the reorder point of -2 tons and 
a lot size of 20 tons provide a locally minimum cost of $18 per 
month . 

In Example 3 .52 we show how to change the basic period. In 
all previous examples the period was one month, or say, 30 days. 

In , Example 3. 52 the period is 6 days. Hence, a demand of 5 tons 
per month is equivalent to 1 ton per day. Similarly, we have a 
carrying cost of $0.2 per ton per 6 days, and a shortage cost of 
$1.8 per ton per 6 days. The replenishment cost of $36 is not changed. 

Note now that the simulation results are identical with the 
long term averages. Also note that all costs are one -fifth of the ^ 
costs in .example 3 * 51 > since the period is now 6 days instead of 
30 days. 




e 



EMANU, CAKHY# SKOHT# KEH.EN 
S# I# 9# 36 



Example 3 . 51 

















assagggy»gg»ig-t' 



uu, 






DEM A . J# CAPrtY# SHO r<T# REPLEN 
? 1 > to. 2# 1*6# 36 



axacple 3*52 



POIN i#LOT SI2E#SIM#INIT# HUN# PHI NTS# AVER# 00 SI S# STEP# DATA 
? -2# 20# 1# A# ltoto# 8# 1# 1# 0. 1# 1 



REORDER POINT = -2 



LOT SIZE = 20 



SIM ULATION 



PER. bEGIN DEMAND END 



1 

2 

3 

A 

5 

6 
7 
0 



A 

3 

2 

1 

0 
- 1 
18 
17 



3 

2 

1 

0 

- 1 
-2 
17 
16 



CARRYING 
3* 5 
2. S 
US 
• 5 
0 
0 

17* 5 
16* S 



SHOHTAGES 

0 

0 

0 

0 

• 5 
1. 5 
0 
0 



REPLENISHMENTS 

0 

0 

0 

0 



1 

0 

0 



FOK A SIMULATION OF 100 PERIODS* 



CARRYING 



SHO RTAGES KEFL ENI SHMENTS 



AVERAGES 
UNIT COSTS 
COSTS PER PERIOD 



8. 1 

• 2 

U 62 



• 1 
1.8 
• 18 



.05 

36 

1.8 



TOTAL COSTS PER PERIOD = 3.6 



THE 



LONG TERM 
CARRYING 



AVERAGES 

SHO RTAGES REPL EN I SHM EN TS 



AVERAGES 
UNIT COSTS 
COSTS PER PERIOD 



8. 1 

• 2 

1.62 



• 1 
1.8 
• 18 



• 05 
36 
1.8 



TOTAL COSTS PER PERIOD * 3.6 



THE TOTAL 



COST 



TABLE 



19.9 



LOT SIZES 
20 



20 . 1 




3. 60065 
3. 60005 
3. 600A5 



3. 0005 
3.6 

3. 6005 



3. 600A5 
3. 6000 5 
3. 60065 



^5 
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The change of period length is also illustrated in 
iScamples 3.53 and 3»5^* The examples deal with an inventory 
system in which the demand is 2^+00 parts per year, the carrying 
cost is $0.56 per part per year, and the replenishing cost is 
$1*2. IIo shortages are allowed. 

One can show that the optimal reorder point is 0 and 
that the optimal lot size is uOO parts. Note the meaningless 
output in the simulation of Example 3*^3 and how this is 
corrected in Example 3»5^» 



DEMAND# CARRY# SHORT# KEPLEN 
? 2400# 0.56# 99999# 42 



Example 3 .$3 



POINT# LOT SIZE# SIM# INI T# RUN# PRINTS# AVER# COSTS# STEP# DATA 
? 0# 600# 1# 600# 10# 5# 1# 1# 10# 1 



REORDER POINT = 0 



LOT SIZE * 600 

M ULArT ION 



PER. 

1 

2 

1 

3 
1 

4 
1 

5 
1 



BEGIN DEMAND END 



CARRYING 



SHORTAGES KEPLEN I SHMENTS 



600 


2400 


-1800 75 


675 


1 


- 1200 


2400 


-3600 


0 


2400 


-3000 


2400 


- 5400 


0 


4200 


-48 00 


2400 


-7200 


0 


6000 


-6600 


2400 


-9000 


0 


7800 



FOk A SIMULATION OK 10 PERIODSi 



iScamplo 3*33 (Cant'd) 



CARRYING 



SHORTAGES REM. EN I SHM EN T S 



AVERAGES 
UNIT COSTS 
COSTS PER PERIOD 



7 . 5 
• 56 
4.2 



8707. 5 

99999 

870741292 



I 

42 

42 



TOTAL COSTS PER PERIOD = *70741338 



THE LONG TERM AVERAGES 



CARRYING 



SHORTAGES REK. EN I SHM EN T S 



averages 

UNIT COSTS 
COSTS PER PERIOD 



300 
• 56 
168 



0 

99999 

0 



4 

42 

168 



TOTAL COSTS PER PERIOD = 336. 



THE 



TOTAL COST TABLE 



590 



LOT SIZES 
600 



610 



R. POINT 
-10 
0 
10 



8804.99 
336* 047 
341. 647 



8663. 7 

336. 

341.6 



8 527. 13 
336. 046 
341. 646 



*************** *4* 

DEM AN D# CARRY # SHO RT # KEPL EN 
? 200# 0. 046667# 99999# 42 



Example 



POINT* LOT SIZE. SIM. INI T. KIN. PRINTS. AVER. COSTS. STEP. DATA 
? 0, 600# 1# 400# 12# 6# t# 1* 10* » 



REORDER POINT = 0 



LOT SIZE * 600 







l 






SIMULATION 



example 3«5^ (Cont'd) 



PER. BEGIN DEMAND 

i 400 am 

a am am 

3 6 m am 

a Am am 

5 200 am 

6 6 m 200 



END 


CARRYING 


200 


300 


0 


100 


400 


S00 


200 


300 


0 


100 


400 


500 



SHORTAGES 

0 

0 

0 

0 

0 

0 



KEPLENI SHMENTS 
0 
1 
0 
0 
1 
0 



FOR A SIMULATION OF 12 



PERIODS! 



CARRYING 



SHO RT AGES RE PL ENI SHMENTS 



AVERAGES 
UNIT COSTS 
COSTS PER PERIOD 



300 

• 046667 
1 A* 0001 



0 

99999 

0 



•333333 

42 

14. 



TOTAL COSTS PER PERIOD « 28.0001 



THE 



I 



LONG 



TERM 



AVERAGES 



CARRYING SHORTAGES REPLENISHMENTS 



AVERAGES 
UNIT COSTS 
COSTS PER PERIOD 



300 

•046667 
14 * 0001 



0 

99999 

0 



• 333333 

42 

14. 



TOTAL COSTS PER PERIOD * 28 •0001 



THE 


TOTAL 


C 0 


ST T A B L 


E 






LOT 


SIZES 






590 




600 


610 


R. POINT 
-10 


8502.03 




8 360. 79 


8224. 18 


0 


28.0041 




28.0001 


28.0039 


10 


28. 4707 




28. 4668 


28. 4706 




* * * * * 4c * * * * * * * * * * * ** * * * 

DEMAND# CARRY# SHORT# REPLEN 
Z STOP 



RAN 



30 SEC. 



lLL 



i 



ERiC 















«^*G 
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h. Construction and Solution of the Model 



A thorough understanding o f how the system behaves is 
necessary before the student can construct a model ol' the system. 
This understanding is achieved by the use- ol tne simulation option 
ol' tne- pi’ogram. 

Once the system is completely understood, one should attempt 
to predict the Ion;' tern averages lor any primary input and any 
reorder points and lot sices. This is the essence ol constructing 
a mathematical model of the system. The student can check hio 
model by comparing numerical resilts of his model with those given 
by t he program. 

Alter the model lias been constructed, the student must find 
an analytical method for the determination or the optimal reorder 
point and lot size. Although he may use tne costs option oi the 
program to test his optimizing technique, the student must 
mat ho : nat i cally justify his results - specific precaution should be 
t alien to eliminate the possibility of assuming that a local minimum 
is the global minimum. 



4.1 Understanding the System 

A good way to understand a system is to see it work. The 
student should therefore select some values for (J,U, and Q and 
graphically describe the progress of the system a;, -i .unction ot 
time, how much inventory is carried (or is in shortage )*. when I 
there a repienislmientV etc. Ills selection ol values should re .l*^ 
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the following situations: 

(a) No shortages are allowed (S 2 0). 

(b) No Inventory Is carried (S + Q £ 0). 

(c) Inventory is carried and shortages are allowed (S < 0 and S ♦ Q > 0) 



Next, he should calculate the three averages: Inventory carried* 
shortages, and replenishments. These should be checked against the 
results that the computer gives in the simulation averages. 



4 .2 Construction of the Model 



The general model of the system may be represented by the 



equation 



C(e,q) = ^ 1^(6, q,r) + Cg Ig(s,q,r) + Cj I^a.q.r) 



(1) 



where 



s «= reorder point 



q = lot size 



C & total cost per unit time 



c^= carrying cost per unit per unit time 



Cg= shortage cost per unit per unit time 



c^s replenishing cost 



I^a average amount carried 



Ig= average shortage 



I^s average number of replenishments per unit time. 



r a demand per unit time 

To find for any s,q, and r, one may wish to consider the 



following cases: 



</f 



% 



, * f 



I 



51 













(a) jjo Cnorta a ro Alloyed 

i 

This case is equivalent to seating timt is in' ini to ly large, 
it snoald dc ouvious Mint s S 0, hence the model in: 

C(s,q) - c 1 I i (3,q,r) + ^(s.q.r) G 25 0 ( Li ) 

The student should now chow that £ is a linear fimetion or s and q, 
and that 1^ is not a function of 3. 

(b) Uo Inventory is Carried 

This case ic equivalent to stating that is infinitely large. 
The model ic now : 

C(s,q) = C 2 I 2 (s,q,r) + c^l^(s,q,r) s + q d 0 (i) 

I should be shown to be a linear function of s and q, and I- to be an 

2 ^ 

in Case (a). 

(c) 1-re scribed Lot Sice 

In this case the lot size is fixed, cay, at « - rt, when t is the 
length of the period. Then 1^ is fixed too. Assume then that c y 0. 
The model becomes 

C(~,q p )-c 1 i i (r.,q p ) + c 2 ('0 













D 






I and Lj should bo dote mined for three distinct situations: 



(1) s *; 0, (2) s -f Op *•» 0, and (3) s -: 0 and s + - 0. The 

results for the first two situations correspond to Cases (a) and 
(b) respectively. In situation (3) one has to show that 1^ and 



I are quadratic functions of s. 
d 



(d) The General Case 

After constructing the models of Cases (a), (b), and (c), 
the general model of .equation (l) can be immediately constructed, 



The model should be stated in the fora: 



C(s,q) 



J 



.a 



^ 0 



s £ 0 s + q - 0 



(5) 



V c u 



+ q 5 0 



-./here C a and C b are as in Aquations (2) and (3) respectively and 



C c is in part as in equation (4). 



As the models are being constructed the student can check 



them with the long term averages option of the program. 'L'he 



examples in the first four illustrations of the previous section 
correspond respectively to Cases (a) to (d). 



H.3 Solution of the Model 

Alter the model of a system has been constructed one generally 
wishes to find the optimal values of the controllable variables. 

The process of finding these values is referred to as the solution 



J e 
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2h 



0 t 4lJ model. Thus we are interested in finding the optimal 
•."a lues ox the reorder point and the lot size. We will reler 



to thxsj values as s q and q Q . 



Xt is a/jain advisable to consider the lour cases described 



in the previous section, 



In Case (a) tne student siioi-Ld show that s^ - 0 and that 



is a function of r, c^, and c'^. 



In Case (c), the student should show that c q is a lunction 



of o , c, , and c Q . 

*p 1 <- 



Finally, in Case (d), he should show that s q is not positive 



md a function of r, c , c Q , and c , and that q Q is another function 



of r, c ± , c 2 , and c^. 



The optimal functions which the student derives may be 
checked us inq the table of total costs of the pro/ 'ram. See the 
details in dxaraples 3»13> 3*22, 3*32, 3»^2» 3«^1» 3»^2, 3*^3» 
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5 . Extensions and Problemc for Solution 



5.1 Extensions 



By changing some oi' the properties ol* the inventory system 



oi' Program ISL-2, a variety of nev systems can be obtained. For 



example, if demand is allowed to vary and follows some given 



probability distribution, we would have the probabilistic reorder 
point - lot size system. This tos been done in Program I3L-3 



which is described in another manual. I/?ss extensive extensions are 



stated in the fora of exercises. 



A. Program I3L-2 assumes tliat shortages are made up - this is the 



classical back-order case. What changes should be made in the 



program for the lost-sales case, in which shortages are not made up 



B. . Set end Program ISL-2 for systems in which leadtime is L periods 
The leadtime should be specified with the primary output. 



C. Change Program I3L-2 so that the controllable variables are the 
scheduling period t and an order level S (instead of the reorder point 
s and the lot sine q). In the new system replenishments occur every 
t periods. The replenishment raises inventory to a level 3. 



D. In Program ISL-2 the lot size q is added to inventory instantaneously 
ilxtend the program to allow the addition to stock to be at a rate p. 



B. ibrtend Program ISL-2 to an inventory system with 2 items. For the 
first item the rate of demand is r^, the carrying cost is c ^ ^ per unit 
per unit time, and the shortage cost is per unit per unit time. 




2 ^ 









J V 



The corresponding parameters i'or the second item are r^, c^, and 
c 22’ ^ ie re P lenistiIfien ‘ t cost is c y Thia cofrt is incurred whenever 
one or both items are replenished. 



5.2 Problems for Solution 

A. In ^cample 3.51 there is a marked difference between the 
simulation results and the long terra averages. 

(1) Explain the reason for the difference. 

(2) What would be the simulation results for initial 
inventories of 15? 16? 17? 18? 19? 

B. Explain the simulation results in Example 3 • 53 • 

C. Find the optimal solution of a deterministic reorder point 
lot size system in which the demand is 6 lb. per month, the 
carrying cost is $0.75 per pound per month, the shortage 
cost is $1.50 per pound per month, and the replenishing 
cost is $5^*00* 

D. Solve Problem C above if the respective quantities are 

25 lb., $9 per pound per month, $l6 per pound per month, and 

$ 288 . 

E. Solve Problem D when the reorder point and the lot size must 
be integer multiples of 20 lb. 




v/'3 



v 















Program ISL-2 



27 



L. ^ADDON AND KlCHAKD SACHF-h* THE JOHNS HOPKINS UNIVERSITY 
THE DEI EKM1N1 STIC KEOKDEk FOINT-LOT SI2E SY STEM" 



j 164 

| 1 6 5 REM PRIMARY INPUT 

| 17o i-rtl.NI T”DEMANU# CAkKY » SHORT* KEPI EN" 

180 INPUT U#C1#C2#C3 
| £20 GO SUB 7 000 

I £30 GOTO 1050 

| 9vv 

1000 kEM SECONDARY INPUT 

10 10 GO 0 U b 7000 

1040 

' 1050 FKIN7"POINT>LOT SI2E# SIM# IN IT# RUN# PRINT S* AVER* CO 51 5# STEP# DATA" 

I 1 U 6 0 1 N P U 7 S#U#Al#Ul#N#Nl#A 2# A 3# J 1 * D 

! ID/DpkInT 

i iuhufkINT 

f 1 09 UFKI N 7 "REO RDER POINT = "S#"LOT SI2E = "D 

f IIJUPkINT 

j 1140PKINT 

f j 1150 IF A 1=1 7rt EN £000 
i 1160 IF A 2*1 THEN 3000 

j 1170 IF A 3=1 THEN 4000 

I 1171 PRINT" **+******♦♦♦♦***♦****♦" 

| 1175 IF 0=1 THEN 165 

| 118 0 GOTO 1000 

| 19*9 

j £000PhIN7" SIMULATION" 

I £UD£ LET A 1=0 

| £004 LET A£=0 

| £006 L ET A3= 0 

j £010 IF N 1 = 0 THEN £0 50 
}. £0£0bh IN 7 

I £0G0PKINT 

[ £040PKIN T "r E k* BEGIN DEMAND END CAhKYING SHOK1AGES KEPLENl* 

FOK 1=1 TO N 
IF Q1< = 0 THEN £059 
IF (01-UX0 THEN 2056 
LET I 1 = L1- U/ 2 
LET 12=0 
GOTO 2061 
LET I ! = Ult2/<2*U) 

LET I 2= C D 1- U> t 2/ ( 2* U) 

GOTO £061 
LET 11 = 0 
LET I £= - L 1 ♦ U/ £ 

LET 0£=L1-U 
IF 0£>S THEN 2066 



| £050 

i £051 

! 2052 

1 £0 53 

I 21 54 

! £055 

J £0 56 

| 2057 

1 £0 58 

j £0 59 

I £060 

k , 2061 

fL_. £062 



[ Li.- 



I 



1 1'0 i\ EM BY 
I 1 U 

1 40rnINT '• 
i 50P.KIN1 
IfcDPXl ^ i 
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ISL-2 CONTINUED 



2063 L E i I 3= I 

2064 LEI t*UH+L i 
206 5 GO TO 2069 

2066 LEI I 3» 0 

2067 LET *=G2 

2069 IF I»N1 THEN 2080 

2070PhiNl I S U 1 J UJ L2# 1 1 # I 2« l 3 

2080 LET Ul=y 

2082 LET AisAH-I I 

208 4 LET A2= A2* I 2 

208 6 LET A3=A3+I3 

2090 NEAT I 

2100 LET I 1 = A 1/N 

2110 LET I 2= A2/N 

2120 LET I 3 s A3/N 

2121 PR INT 

2122 PRINT 

2123 PRINT 

2125 PHI NT " FOR A SIMULATION OF "Nf •• PERIODS*" 

2130 GO SUB 6000 
2320 GOSUB 7000 
2330 GOTO 1160 
2999 

3000PKIN T M THE LONG TERM AVERAGES" 

3002 PRINT 

3003 GOSUB 5000 
3005 GOSUB 6000 
3100 GOTO 1170 

3999 

4000 RES'* COST TAbLE 
4002 LET S9*S 

4004 LET 09= 0 
4006 LET SC 1 > = S9 - J 1 
4010 LET SC 2) = S9 
4020 LET SC 3) = S9+ J 1 
4030 L ETGiC 1 ) = 09 - J 1 
4040 LET QC 2) = 09 
4050 LET 0C3) = G9+J1 
4060 GOSUB 7004!) 

4070 PRINT" THE TOTAL COST TABLE" 

4080PRINT 

4090PRINT 

41A0PRIN7" LOT SIZES" 

4T50PRINT" ", G< 1>#Q(2>,G<3> 

4160PKINT"R. POINT" 

4170 FOR 1=1 TO 3 
4130 LET S= SC I ) 

4190 PRINT S» 

4200 FOR J= 1 TO 3 
4210 LET G=U(J> 
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I5L-2 CONTINUED 
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4220 GO SUB 5000 

4230 Ph INI Cl*Il + C2*I2*(;3*i3# 

424® NEAT J 
42 b® PKINT 

4260 NEXT I 

4261 PKINT 
4270 GO TO 1 171 

4999 

5000 KEM I 1# I 2# I 3 
500S LET I 3=U/Q 

5010 IF S< 0 THEN 5050 
5020 LET I l=S+U/2 
5030 LET 1 2=0 
5040 KETUKN 

5050 IF S<-U THEN 509 0 
5060 LET I 1 = < U+S> t 2/< 2*Q> 

50T0 LET I2=St2/C2*U> 

508 3 KETUKN 
5090 LET I 1 = 0 
5100 LET I 2=- S+U/2 
5110 KETUKN 

5999 

6000 KEM AVEKAGES AND COSTS 
6005 GO SUB 7000 

6010PKINT** CAKKHNG SHOKTAGES 

6020 PKINT 

6040PKINT ,, AVEKAGES"# ” *'# I l# I 2# I 3 

6050PKIN T "UN IT COSTS**# ** **# Cl# C2# C3 

6060PKINT"COST S PEK PEKIOD**# C1*I 1# C2*I 2# C3*I 3 
6070PKINT 

6080PKIN T*'TO TAL COSTS PEK PEKIOD = *'C1*I 1 + C2*I 2+C3*I 3 

6170 KETUKN 

6999 

7000KEM 3 SPACE SUB 

7010 PKINT 

7020 PKINT 

7030 PKINT 

7040 KETUKN 

7999 

9999 END 
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I. Introduction 
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This laboratory deals with System 3 which is an extension of 
System 2. Whereas in System 2 demands are known and constant, in 
System 3 demands are probabilistic. Otherwise the two systems are 
identical. As in ISL-2, the program can be used as follows: 

A. For any reorder point and lot size the student can 
find out how the system behaves over a reasonable 
number of periods, and what the corresponding averages 
are. This is a way to determine whether he understands 
the system. 

B. The student can next attempt to build a model for the 
expected behavior of the system in the long run. The 
The construction of the model will probably prove to be 

a challenge. The student can use the computer program to 
cest his model by comparing numerical results that he 
derives on paper with the numerical results derived on the 

computer . 

C. Once the model of the system has been determined, the 
student addresses himself to the problem of optimization: 




How to find the optimal values of the reorder point and the 
lot size . This is strictly a mathematical problem for which 
it is rather difficult to get a closed form solution. The 
student must attempt to find an algorithm leading to the 
desired values. Any numerical results that he obtains can be 
checked with the third type of output available from the 
the computer program. 

-i ■ . v , . 








2. The Probabilistic Reorder Point - Lot Size System 



2.1 Demand s 

The demand X during a reviewing period is a random variable 
which can take discrete values; 0, U, 2U, etc., up to the maximum 
demand of X2. The probability of demand, P(x), must be given. The 
student supplies the value of U, X2, P(0), P(U), . . .,P(X2). 

2.2 Unit Costs 

The cost of carrying inventory is Cl per unit per reviewing period. 
The cost of shortage is C2 per unit per reviewing period. The cost of 
replenishing is C3 per replenishment. The student supplies the values 
of Cl, C2, and C3. 

2.3 Output 

The program allows for three types of outputs : 

A. Simulation . For a run of size M the inventory system is 
simulated by giving M demands from the probability 
distribution supplied by the student. The simulation provides 
detailed printouts for the first T periods. For each of 
these periods the average amount carried, the average shortage, 
and an indication of whether there was a replenishment are 
supplied. At the end of the simulation the overall averages 
and the corresponding costs are printed. The student has to 
supply the values of the initial invertory 10, the run size M, 
and the number of desired printouts T. 



B, Long Term Averages. If this option is indicated the 



computer prints the expected values of the average 
inventory, average shortage, average number of replenish- 
ments per unit time, and the corresponding expected costs. 

C. A Table of E^cted Total Costs . For this option expected 
costs in the neighborhood of the specified reorder point 
and lot size are printed. 

2.4 Decisions 

The student can control the reorder point SI and the lot size Q. 

These variables, as well as the variable s controlling the output, are 
supplied by the student as input during the running of the program. 

When the third output option is elected the expected total costs are 
given for the values of Sl-U, SI, and Sl+U, and Q-U, Q, and Q*HJ. 

2.5 Objectives 

The student supplies the unit costs and the probability distribution 
of demand in data statements. Thereafter, he has three major objectives: 

A. To check whether he understands how the system behaves using 
the simulation option of the program. 

B. To construct the model of the system and to use the model 
to determine numerical values of the expected average 
inventory which will be carried, the expected shortages, the 
expected number of replenishments, and the expected total 
cost of the system. The long term averages option of the 
program allows him to check whether his model indeed 
describes the system. 

Lr o 




C. To develop an algorithm for finding the optimal values 



| of the reorder point and the lot size . The table of 

j expected total costs can be used to check the algorithm. 



3. Presentation of Data and Inputs 
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3.1 Unit Costs: Cl, C2, C3 

Three unit costs have to be supplied by the student. The unit 
carrying cost is in dollars per unit per reviewing period. The unit 
shortage cost is also in dollars per unit per reviewing period. The 
unit replenishing cost is in dollars. The general data statement is: 

9010 IMA < Cl, C2, C3 > 

For example, if the unit carrying cost is $5 per ton per week, the 
unit shortage cost is $50 per ton per week, and the unit replenishing 
cost is $40, the data statement reads: 

9010 TATA 5, 50, 40 

3.2 Probability Distribution of Demand: P(0), P(U), . . . ,P(X2) 

Demand can take the values of 0,U,2U, . . . ,X2-U,X2. The corresponding 
probabilities are ?(0), . . . ,P(X2) . The general data statement is in form: 
9020 I&TA< U,X2,P(0),P(U), ...,P(X2) > 

For example, for demands of 0,2,4 ,6, and 8, with corresponding probabilities 
of .05, .24, .38, .21, and .12, the data statement is: 

9020 I&TA 2,8, .05, .24, .38, .21, .12 




Obviously, the sum of the probabilities must be equal to 1. The cumulative 
probability distribution is computed by the program and is printed 
immediately. If the cumulative probability distribution of X2 is not equal 
to 1, the student should stop the execution of the program and check his 



data statements, 
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3*3 Inputs 

The program requires 8 inputs during its running: the reorder point 
SI, the lot size Q, an index Y related to simulation, the initial 
inventory 10, the simulation run size M, the number of periods to be 
printed T, an index S related to long term averages, and an index R 
related to the table of expected costs. 

Each of the indices Y,S, or R should be either 0 or 1. 

When Y=0, the simulation option is not used. In that case, the 
values of 10, M, and I are irrelevant; however, they must be supplied. 

If Y=l, then the simulation starts with an initial inventory of 
TO units, it is carried out over M periods, and the details of the 
first T periods are printed. If T=0, only the averages for the 
simulation are printed. 

If S=l, then the long term averages are printed. (These can be 
compared with the simulation averages.) If S=0, then the long term 
averages are not printed. 

If R=l, then a table of expected costs is printed in the neighbor- 
hood of the reorder point SI and the lot size Q. Otherwise such a table 
is not printed. 

The request of inputs is in the forrn^ 

POINT, LOT SIZE, SIM, EOT, RUN, PRINTS, AVER, COSTS 
? < SI , Q , Y , 10 , M , T , S , R > 



For example, the following inputs may be supplied: 

POINT, LOT SIZE, SIM, INIT, RUN, PRINTS, AVER, COSTS 

? 0 , 10 , 1 , 10 , 100, 5 > 1 > 1 




(> V 
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In this example the student is interested in a reorder point of 0, and 
a lot size of 10. He wants a simulation with an initial investing of 
10 to be run over 100 periods. He asks that the details for the first 
5 periods be printed. He also wants the long term averages and a table 
of expected costs. 
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4. Illustrations 

Two illustrations are given in this section. The first illustrates 
the. full capabilities of Program ISL-3* The second deals with a 
deterministic reorder point - lot size system and allows a comparison with 
Program ISL-2. 

4.1 Illustration 1 

The reviewing period in a reorder point - lot size system is one week. 

The carrying cost is $5 per pound per week, the shortage cost is §50 per 
pound per week, and the ordering cost is $40 per order. 

Demand may be assumed to occur in units of 2 pounds, with a maximum of 8 
pounds in any week. The demand probabilities may be assumed to be P(0) = .05, 
p(2) = .24, P(4) = .38, P (6) = .21 and P(8) = .12. 

For some reorder points and lot sizes one wishes to study the behavior 
of the system and to check long term averages. One also wants to check whether 
some specified reorder point and lot size give a minimum total expected cost. 

The data, the computer outputs, the inputs, and the results are given 
on Pages 8 to 10. 
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9U1U DATA 5#5J#4J 
9J2U OATrt 2,3, .J5# 
RUM 



I# • 0*5 P • 



# • 



I SL - 3 * 



21836 W2 WED 10/25/67 



THE PROBABILISTIC REORDER PO I NT- LOT SI2E S/STEM 



CARR/ 1 MG COST = 5 PER UMIT PER PERIOD 

SHORTAGE COST * 50 PER UMIT PER PERIOD 

REPLEM I SR I MG COST = 40 PER REPLEM l SHMEM T 



DEM AMD 
0 
2 
4 
6 
3 



PROBABILI T/ 

• 05 

• 24 

• 33 

• 21 
• 12 



CUM UL ATI VE 

• 05 
.29 

• 67 
.83 
1. 



♦ *♦**^***4: **:•:****: ft*******:*#* 



PO I M T# LO TS 1 2 E# SIM# IMIT# RUM# PRI M TS# AVER# COSTS 
? -4# 14# 1# 10# 10# 10# 0, 0 



Example h.ll 



REORDER POIMT =- 4 



LOT SI2E = 14 



SIM'JLATIOM 



PER 8EGIM DEM AMD EMD 



CARR/ l MG 



SHORTAGES REPLEM I SHMEM TS 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 



10 

4 

2 

0 

10 

8 

4 

2 

3 

4 



6 

2 

2 

4 

2 

4 

2 

3 

4 
2 



4 

2 

0 

-4 

8 

4 

2 

-6 

4 

2 



7 

3 

1 

0 

9 

6 

3 

• 25 
6 
3 



0 

0 

0 

2 

0 

0 

0 

2.25 

0 

0 



0 

' 0 
0 

1 

0 

0 

0 

1 

0 

0 




FOR A SIMUL ATIOM OF 10 



PERIODS: 



DEM AMD 



CARR/IMG 



SHORTAGES REPLEM I SHMEM TS 



AVERAGES 
UMIT COSTS 
COSTS PER PERIOD 



3.6 



3.825 

5 

19. 125 



• 425 
50 

21.25 



• 2 
40 
8 * 



TOTAL COST PER PERIOD = 43.375 



******** ****$**:>:*& 4c ******:(: $* it**:*******:*^ 4c,* ******** * ********** 
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P0INT*L0TSI3E* SI ^ ^ I M I RUN* PRI M TS# AVER* COSTS 

? -4* 14* 1* 10* 500* 10* 1* 1 



REORDER POINT =-4 



LOT SI2E s 14 



SIMULATION 



PER BEGIN DEMAND END 



CARRYING 



1 

2 

3 

A 

5 

6 

7 

8 

9 

10 



10 

8 

2 

•2 

6 

2 

8 

0 

8 

4 



2 

6 

4 

6 

4 

8 

8 

6 

4 

4 



8 

2 

-2 

-8 

2 

-6 

0 

-6 

4 

0 



9 

5 

• 5 

0 

4 

.25 

4 

0 

6 
2 



FOR A SIMULATION OF 500 PERIODS: 



DEMAND 



CARRYING 



AVERAGES 
UNIT COSTS 
COSTS PER PERIOD 



4. 172 



2.38067 

5 

1 4. 4033 



TOTAL COST PER PERIOD = 76.4567 



LONG TERM AVERAGES 



DEMAND 



CARRYING 



AVERAGES 
UNIT COSTS 
COSTS PER PERIOD 



4.22 



2.91571 

5 

1 4. 5736 



TOTAL COST PER PERIOD = 77.9214 



LOT SI2E » 
R. POINT 
-6 
-4 
-2 



THE EXPECTED TOTAL COST TABLE 
12 14 



130. 342 
84. 3333 
56. 6533 



115.929 
77.9214 
55. 6236 



-LIT 



Example 4.1 2 



SHORTAGES REPLENI SHMENT5 



0 

0 

. 5 
5 
0 

2.25 

0 

3 

0 

0 



1.00267 

50 

50. 1333 



1.02571 

50 

51.2357 



16 



106.369 
74. 3625 
56. 1062 



0 

0 

0 

1 

0 

1 

0 

1 

0 

0 



SHORTAGES REPLENI SHMENTS 



.293 

40 

11.92 



SHO RTAGES REPL EN I SHM EN TS 



. 301429 
40 

12.0571 



* * ***********************.:* ********************************** * 

^ n t - T ft ; — — — 

■ERIC 




i 














I 



I 



1U 



POI 'JT#LOTSl?E# SIM, INIT# RUN# PRI N TS# AVER* COSTS 
? £, 3, 0# 0* O# 0* U# 1 



REORDER POINT = 2 



LOT SI n E = 8 



J^ample 4.13 



THE EXPECTED TOTAL COST TABLE 



LOT SI 3 E = 
R. POINT 
0 
2 
4 



6 


3 


10 


53. 5833 


43. 75 


46. 89 


50.0167 


48. 575 


48. 75 


56. 5333 


55.9625 


56.66 



****** ********************** ********************************* 

Example 4«l4 



POINT# L0TSI2E# SIM# INI T# RUM# PRINTS# AVER# COSTS 
? 0, 10, 1# 4# 1000# 0# 1# 1 



REORDER PO I N T = 0 



LOT SI2E = 10 



FOR A SIMULATION OF 1000 

DEMAND 



SIMULATION 

PERIODS: 
CARRYING 



SHORTAGES REPLENI SHMENTS 



AVERAGES 
UNIT COSTS 
COSTS PER PERIOD 



4. 336 



3.97667 

5 

19.8333 



.214667 

50 

10. 7333 



. 434 
40 

1 7. 36 



TOTAL COST PER PERIOD = 47.9767 

LONG TERM AVERAGES 

CARRYING 



averages 

UNIT COSTS 
COSTS PER PERIOD 



DEMAND 
4. 22 



SHORTAGES REPLENI SHMENTS 



4.032 

5 

20. 41 



• 192 
50 
9. 6 



• 422 
40 

16.38 



TOTAL COST PEP PERIOD = 46.89 

THE EXPECTED TOTAL COST TABLE 




3 


10 


12 


67. 7625 


60. 1 


56.6533 


48. 75 


46.89 


47.3167 . 


43. 575 


48. 75 


50. 5333 



[ERIC 
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Exaiaples 4.11, 4.12, 4.13, and 4.l4 illustrate the following 

points : 

a. The simulated demands follow the given distribution. Different 
simulations use different sequences of demands (see Examples 4.11 
and 4.12). The mean demand is 4.22 (see Example 4.12 under 
Long Term Averages). However, the simulated mean demand varies 
even for large simulations (in Example 4.12, for a simulation of 500, 
it is 4.172, while in Example 4.l4, for a simulation of 1000, it is 
4.336). 

b. The method of obtaining simulation averages and costs is illustrated 
in Example 4.11, where only 10 periods are simulated. This example 
also clearly illustrates the behavior of the system from period to 
period, the resulting average amounts that are carried and/or are 
short, and whether replenishments are ordered (and delivered). 

c. The method for computing long term averages is not illustrated, of 
course. This is one of the principal tasks the student has to work on. 

d. Examples 4.12 and 4.l4 can be used to compare simulation averages and 
long term averages. 

e. Example 4.13 provides an illustration of a reorder point SI and a 
lot size Q which are, in a certain sense, locally optimal, since 



C(S1,Q) £ C(S1 t u,Q) and C(S1,Q) £ C(S1,Q ± U) 



As can be seen from Example 4.l4, the optimal solution appears to be 
S1=0 and Q = 10, with a minimum cost of $46.89 per week. 




6 -? 












4.2 Illustration 2 



Program ISL-3 can handle both probabilistic and deterministic 
reorder point - lot size systems. However, as compared with program 
ISL-2, it is not flexible for deterministic systems in the display of 
The Expected Total Cost Table. In Program ISL-2 the step for various 
reorder points and lot sizes can be specified with the input . In 
Program ISL-3 the step equals the basic demand unit U. 

Examples 4.21 and 4.22 illustrate a system with deterministic 
demand of, say, 5 tons per month. The expected costs in Example 4.21 
should then be interpreted to be in dollars per month. If one desired 
a s malle r basic quantity unit, say, 1 ton, then the basic time unit 
would have to be 30/5 * 6 days. For this time unit the data and 
results are as in Example 4.22. The costs are now in dollars per 
6 days. Thus the cost for a reorder point of -2 and a lot size of 
20 is $3.6 x 5 = $18 per month. 

Construction and Solution of the Model 
The properties of the inventory system to be studied can be 
summarized as follows: Demand x during some reviewing period w is 
uniform and follows a discrete probability distribution P(0), P(u), 
p(2u), •••» 1> ( x taaac ) where u is some basic unit and is the maximum 
demand during w. The carrying cost is c^ per unit per unit time, the 
shortage cost is Cg per unit per unit time, and the replenishing cost 
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9010 DATA 1 >9,36 
9020 DATA 5 >5,0,1 
HUN 



Example 4.21 



ISL-3* 21:55 W2 WED 10/25/67 

THE PROBABILISTIC REORDER POINTt LOT SIZE SYSTEM 



CARRYING COST = 1 PER UNIT PER PERIOD 

SHORTAGE COST = 9 PER UNIT PER PERIOD 

REPLENISHING COST = 36 PER REPLEN I SHMENT 



DEMAND 

0 

5 



PROBABILITY 

0 

1 



CUMULATIVE 

0 

1 



5 }c jJc rj: ^ i}; ;{; ^ ?}c ifc sjc if: ^ ^ % s}c A sfc sfc ^ sjs ^ ^ '*£ ^ ^ ^ ^ *5* ^ *5* ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ **»* ^ *V 



POINT, LOTS I ZE , S I M , I N I T , RUN , PR I NTS > A VER > C 3ST S 

? 0 > 20 , 0 , 0 , 0 , 0 , 0 , 1 



REORDER POINT = 0 



LOT SIZE = 20 



LOT SIZE = 


15 


20 


25 


R. POINT 
-5 


22.8333 


20.25 


19.7 


0 


1 9.5 


19 


19.7 


5 


24.5 


24 


24.7 



****************************************** **************** 




THE EXPECTED TOTAL COST TABLE 



P0 1 NT, LOTS I ZE, SIM, INI T,RUN, PRINTS, AVER, COSTS 
?• STOP 



6 SEC. 
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Example 4.22 



9*51 a DATA .2* l .8* 36 
9 ') ‘20 DA i A 1 # 1»! : :»1 



ISL-3* 22:00 >'.2. .*.‘E0 13/25/67 

THE PROBABILISTIC r£OrDE:\ POINT- LOT 51 ?£ SYS TEX 



CA r\ ri V 1 N G COj».I = • 2 

Sri Oh I A fit: COST = 1.8 
HEPLEN 1 5H I. rifi GOST “ 36 



P £.-< UiV IT P £r P£*l j 0 
PEr UN' I T rEn PErIOm 
PEH P.EPLEN I SHX'Ei\ T 



BEX ANU 
0 

1 



PkOBABIlI 1 y 
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1 



CU/j JLA fWi 
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l 



£ 5$; Jr £ |c f; :{t fr.y 4. 4; .Js *,: J: f: * :f« sjf if* *j« * •• ?■ 4 = ¥ * ¥ S* ¥ S ¥ *S ❖ -t- r : * 1: *= $ +' ^ *+' +' ■> y * * 



P<> I i'i T p L. «..) I r I r. E * S I <-> > 1 1 \ 7. I * P u v * P -v i .V IS# A UEs# GO S i S 



? -2 






, / ’j f* 
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LOT SI EE = 20 



THE EXPECTED TOTAL COST TABj-E 



LOT SIS: 

r. POINT. 

-3 

-2 

-1 



19 



20 



2t 



3. 668 42 
3.60526 
3.64737 



3.65 
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3*65 



3.64286 
3 • 60 47 6 
3.6619 
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An(s,q) policy is used. That is, inventories are reviewed 
every period w. When the amount on hand is s or below, a replenish- 
ment of q. units is added to inventory immediately. If the total 
amou nt is st ill at s or below, another q units are added. As 
many replenishments of q are added until the amount on hand is 

! larger than s. However, only one replenishment cost of c~ is incurred, 

| • 3 

| no matter how many q's are needed to raise inventory above s. 

f 

l The model of the system can be represented by: 





$ 



I 



C(s,q) = c^I^(s,q) + c 2 I 2 (s,q) ^ 

where C is the expected total cost of the system per unit of time, 

1^ and Ig are respectively the average amount carried and the average 
shortage, and 1^ is the number of replenishments per unit time. 

The solution of the model is a pair of optimal values s q , q Q 
which minimise C. Namely 



0(8^^) S C(s,q) 

5.1 Construction of the Model 

To find C(s,q) one only needs to find I^(s,q), I 2 (s,q), and I^(s,q). 

It is suggested that the student first find I^fsjq) . He should then find 
J-(s,q) when I 2 (s,a) = 0. Next he should find I^Sjq) - I 2 (s,q). Finally, 
he should find Ig(s,q) . 
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As the student proceeds in finding the appropriate part of the 
model, he should use the Long Term Averages option of Program ISL-3 
to check his results. To do so he'll have to assume numerical values 
for w, u, x jn3x , P(x) (x=0, . ..jX^^), c 1> c^, and c^. He'll then have 
to prepare the corresponding parameters for the data in Lines 9010 and 
9020. He should next compute 1^, Ig, and I^> for several sets of 
suitable values of s and q. He will then be ready to use ISL-3 to 
check whether his answers agree with those provided by the program. 

The student may also wish to use the Simulation option of the 
program to help in understanding the behavior of the system. For 
example, he may wish to examine the probability distribution of the 
amounts on hand at the beginning of each reviewing period. 




This manual contains the detailed listing of Program ISL-3. 
Naturally, the model of the system can be inferred from this listing. 

It is hoped that the student will refrain from doing so. He should 
attempt to build the model only through his understanding of the 
inventory system. Program ISL-3 should only be used to check numerical 
results, in the same manner that a physicist checks a model of a 
physical system when he performs an experiment. 
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5.2 The Solution of the Model 

The problem of finding the values of s and q which minimize the 
expected total cost of the system C(s,q) is an optimization problem. 

The student is required to develop an algorithm which will yield the 
optimal values s q and q Q . This algorithm is not part of Program ISL-3 • 
The only thing the program can do is display the costs in the neighbor- 
hood of some specified s and q,as in Examples 4.12, 4.13, and 4.l4. 

It should be pointed out that the function C(s,q) should not be 
assumed to be convex. Special precaution must be taken to eliminate 
the possibility of assuming that a local minimum is the global minimum 
(e .g. , compare the expected total cost tables of Examples 4.13 and 4.l4). 

Insights int o the optimization problem may be gained by studying, 
in order. Models A to E: 

Model A. The units of shortage cost and replenishing cost are 
relatively very large compared to the unit carrying cost. 

Model B. The unit of shortage cost is relatively large compared 
to the units of carrying cost and replenishing cost. 

Model C. The lot size is fixed. An algorithm is required to 
find only the optimal reorder point s q . 

Model D. An algorithm is required for finding a local minimum 
of C(s,q) . 

Model E. An algorithm is required for finding the global minimum 
of C(s,q). 

7 ^ 
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6. Extensions and Problems for Solution 



6.1 Extensions 

Program ISL-3 may be extended by relaxing or changing some of the 
properties of the inventory system: 

A. Instead of zero leadtime one may have a leadtime of L periods. L may 
be a constant or a variable with a probability distribution G(l). 

B. The backorder assumption is replaced by the lost sales assumption. 

That is, shortages are not made up and instead of Cg being in dollars 
per unit quantity per unit time, it is in dollars per unit quantity. 

C. The reviewing period w is a variable subject to control. 

D. A replenishing cost of c^ is incurred for each lot size (instead of 
one cost for each replenishment - see Page 15 ) . 

6.2 Problems for Solution 

A, man inventory system with an (s,q) policy, no shortages are allowed. 
The amounts in inventory are reviewed every 2 weeks. The probability 
distribution of demand during a 2 -week period is given by P(0) = 0.25 j 
P(5) = 0.20, P(10) = 0.10, P(15) = 0.20, and P(20) = 0.25. The 
carrying cost is $3,20 per unit per week. The replenishing cost is 

$180.00. 

Assuming that the lot size must be a multiple ox 5 units, 
find the optimal lot size and the corresponding expected minimum 
total cost of the system. 
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B. In a probabilistic lot size system with no shortages, the 
probability density of demand during a reviewing period w is 
f(x) = 6x(l-x), 0 ^ x £ 1* 

Find the optimal lot size as a function of A = c^/c^w. 

C. in a reorder- point system the prescribed lot size is 6 units. 

The reviewing period is 1 week. The probability distribution 

of demand during the reviewing period is P(0) = 0.05, P( 2 ) = 0.2k, 
P(lf) = 0.38, P(6) = 0.21, and P(8) = 0.12. The unit carrying 
cost is c x = $5 per week. Find the optimal reorder 
point for three possible unit costs of shortage Cg = $5> Cg = $50, 
and Cg = $100 per week . 

D. The probability density of demand during the reviewing period 
is f (x) = 6x(l - x) 0 £ x £ 1. An (s,q) policy is employed. 

The numerical values of the units of carrying cost and shortage 
cost are equal. What are the optimal reorder points for lot sizes 

of 0.1 and 2.0? 

E. Is the solution found in Example U.l^ the optimal solution? 

F. Solve an extension of Problem D when the reviewing period is 

2 weeks, the carrying cost is $5 per unit per week, the shortage 
cost is $5 per unit per week, and the replenishing cost is $20. 

2 

G. Solve Problem B when f(x) = I2(x - l/2) , 0 £ x £ !• 




H. Find the solution of the probabilistic reorder point - lot size 
system when w = 1, P(0) = 0.08» P(lO) = 0.10, P(20) = 0.20, 
p(30) a 0.30, P(40) a 0.16, p(50) a 0.10, and p(60) = 0.06, 

o^a 1, Cg = 10, and = 25. Show that the minimum expected cost 
of the system is 49*35* 

I. Solve an extension of Problem A when the unit cost of shortage 



is $50 per we*k. 



300OREM BY E.MADDOR AMD I.PRESSMAM 

3010 DIM PC 50>#FC50>*GC50># GC 50># WC50>#MC 50)# VC 50) 

3020MATP=3ERC 50) 

3030MATQ=2ERC 50> 

3040 PRIMT” THE PR08ABILI STIC REORDER POIMT- LOT SI2E SYSTEM” 

3050 PRIMT 
3063READCWC2# C3 

3070PRI MT”CARRY IMG COST =”C1S”PER UMIT PER PERIOD” 

3030PRIMT”SHORTAGE COST = ”C2S ,,p ER UMIT PER PERIOD” 

3090PRIMT”REPLEMI SHI MG COST s”C3l”PEK REPLEMI SHMEMT” 

31 00LETW1 = 1 

3110 PRIMT 

3120 READ Up <2 

3130 LET J2=IMTCX2/U+.0i > 

3 1 35LETX 5=0 

3140 FOR J =0 TO J2 

3150 READ PCJ) 

3155LETX5=X5+J*U*PCJ> 

3160 MEXT J * 

3170 LET .J* J2 

3130 LET A=0 

3190 LET 3=PCJ)/J 

3200 LET QCJ)=A+3/2 

3210 LET A=A+3 

3220 LET J=J-1 

3230 IF J=0 THEM 3250 

3240 GOTO 3190 

3250 LET QC0)=PC0)+A/2 

3260PRI M T” DEM AMD'S SS PR0BA3ILI TY”> "CUMULATIVE” 

3273 LET VC0> = PC0) 

3230 LET GCU>=QC0> 

3290 LET WC0>=GC0) 

3300 LET FC 0 ) =PC 0) 

3310 LET M C 0 > - UC 0) 

3323F0R J=0TQ 49 

3330 LET FCJ+1 )=FC J)+PCJ+1 ) ! 

3340 LET VC J+ 1 )=VC J)+FC J+l ) 

3350 LET GC J+l > = GC J)+QCJ+1 ) 

3360 LET WCJ+1)= VKJ> + GCJ+1> 

33651FJ> J2THEM3333 
3370PRI MTJ^U> PCJ>#FCJ> 

3380 LETMC J+1)=MCJ> + WCJ+1> 

3390 MEXT J 
4003 PRIMT 

40 10 PRI M T”#** ******* * < s * , * c: f £ f 8 4« 4s**-#**** is* * fe ^ 4s ♦ -?e ♦ ♦ * ****** -f: tJ: - fr -Jc >!: It i: 

4020 PRIMT 

4030PRI MT”PO I M T»LO TSI SIM> IMIT> RUM^PRIMTS* AVER# COSTS” 

4043IMPUT SWQ»Y*I0*.M#T#S#R 
4045PRIMT 

4046PRIMT -n-i 
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4047PRINT"REORDER POINT ="Sl*"LOT SIcE = "Q 

4043PRINT 

4049PRINT 

4050IFY= 1THEN5000 

4060 1 FS= 1 THEN 6 0O0 

4070 1 F. H= 1 THEN 7030 

403 060 TO 4000 

5000 PRINT 

5010 PRINT 

5020 PRINT" SIMULATION" 

5030 PRINT 
5030IFT< = 0THEN5120 
5090 PRINT 

5100 PRINT"PER BEGIN DEMAND END "# "CARRYING"* "SHORTAGES 

5110 PRINT 

5120 LET 91=10 

5130 LET L 1=0 

5140 LET L2=0 

5150 LET L3=0 

51 55LETL6=0 

5160 LET K=0 

5170 FOR J= 1 TO M 

5130 LET Y= RDMC 0> 

5190 LETJ9=0 

5200 IF FC J9 ) > Y THEN 5230 

5210 LET J9= J9+ 1 

5220 GOTO 5200 

5230 LET X = J9*U 

5240 LET El= Bl-X 

5250 IF 8 1 < = 3 THEN 5300 

5260 IF El <0 THEN 5330 

5270 LET II =C31+El>/2 

5230 LET 12=0 

5290 GOTO 5350 

5300 LET 12= C-E1-BD/2 

5310 LET I 1 = 0 

5320GO TO 5353 

5333 LET I 1 = B1 1 2/ C 2* C 81 - El > > 

5340 LET I2=E1»2/C2*C31-E1>> 

5350 IF El > SI THEN 5410 
5360LETI 3= 1 
5370 LET E2= El 

5380 LET E2=E2+Q . 

5390 IF E2> 31 THEN 5430 

5400 GOTO 5330 

5410 LETI 3=0 

5420 LET E2= El 

5430 IF J> T THEN 1 450 

5440 PRINT JJ BUX* El# II# 12# 13 

5450 LET L 1 =L 1 + 1 1 
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ISL-3* CONTINUED 



! 5460 LET L2=L2+ 1 2 

|T 5470 LET L3=L3+ 1 3 

5475LETL6=L6+X 
5480 LET B1 = E2 
5490 NEXT J 
5500 PRINT 

5510 PRINT “FOR A SIMULATION OF"MS" PERIODS:" 
5520LETI 1=L1/M 
5530LETI 2=L2/M 
5540LETI 3=L3/M 
5545LETX6=L6/M 
5550GO SUBS 500 
5560GO TO 4060 
6000GOSU33000 
6010 PRINT 
6020 PRINT 

6030 PRINT" LONG TERM AVERAGES" 

6035LETX6=X5< 

6040GOSU33 500 
6070GO TO 4070 
7000 PRINT 
7010 PRINT 

7020 PRINT" * THE EXPECTED TOTAL COST TABLE" 

7030 PRINT 

7040PRINT" LOT SI2E = "# Q-U# Q* Q+U 
7050PRINT"R. POINT" 

7060 LET S1 = S1-2*U 
7070 LET Q=Q+U 
7030 FOR 1=1 TO 3 
7090 LET S1 = S1 + U 
7100 PRINT Sl> 

7110 LET Q=Q-3*U 
7120 FOR J=1 TO 3 
7130 LET Q=Q+U 
7140 GOSUB 8000 
7150 PRINT C, 

7160 NEXT J 
7170 PRINT 
7180 NEXT I 
7 390GO TO 4000 

8000 LET B= INTC C Sl + Q-U)/U+* 01 > 

80 10IFB>=0THEN3060 
8020LETI 1=0 
8030GO T08 1 50 

8060 LET A= INT((S1-U)/U+.01) 

8070 IF A<0 THEN. 8100 
8030 LET I l=C<Ut2)/Q)*CN<B>-NCA>> , 

8090GO T03 1 50 

8100 LET Il = ((Ut2>/Q)*N(B) 

3150 LET 12=1 l+X5/2-((Q+U)/2>-Sl . 
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ISL-3* COMTIMUED 



6152 If 1 2> l E- 6 T4EM 8160 
8154 LET 12=0 

6160 LET D= INT<CQ-U>/U+.01> 

8200 LET I3=C 1-CU/Q>*VCD> J/Wl 
8210 LET C=C1*I1+C2*12+C3*I3 
8220RETURNJ 
8 500PRIMT 
8 502PRIMT 

5 505PRIMT" "# "DEMAND’S "CARRTI'IG"# "SHORTAGES REPLEMI SH'IEMTS" 
8510PRINJT 

8520PRINT"AVERAGES M # X6# II# 12# I 3 

8 530?RI MT"UNI T COSTS'*# '* "#C1#C2#C3 

8 540PRIN T"CO STS PER PERIOD "# Cl *1 1 # C2*I2#.C3*I 3. . 

8 5S0PRIMT 

8 560PRt > J T" TOTAL COST PER PERIOD = "C1 *1 1 + C2*I 2*C3*I 3 
8 570RETURM 

9000 RE^ DATA <C 1 # C2# C3# U# X2# PC 0 ) # PC U) # • • • PC X2) > 

9310 DATA 5# 50# 40 

• 9020 DATA 2# 8# *05# *24# • 33 # • 2 1 # • 1 2 
9999EMD ~ ' 
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1. Introduction 



This laboratory has been designed for the analysis of a variety of 
inventory systems* In particular it is suitable for 



A* Construction of mathematical models. 



B. Checking of Decision Rules. 

C. Study of sensitivity of costs to parameters and/or decision rules. 

D. Comparing methods of generation of random demands and their use. 

i 

3. Analysis of complex inventory systems for which useful mathematical models 
cannot be constructed. 



2. Properties and Policies of the Available Systems 



2.1 Demand 



Demand during a reviewing period may be deterministic or probabilistic. 
It occurs uniformily over the period. The characteristics of the demand are 



supplied as data. 



2.2 Replenishments 

Replenishments occur only at the beginning of reviewing periods. 
Quantities ordered are always delivered after the lapse of the appropriate 



leadtime . 



Leadtime may be deterministic or probabilistic. Its characteristics are 



supplied as data. 
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2.3 Costs 

Four types of costs can be balanced: carrying, shortage, replenishing, 
and reviewing. Information needed to determine the unit costs is supplied 
as data. 

Carryi ng costs are based on the unit cost of the item and on the 
percentage of annual carrying cost. The unit is in [$]/[Q][T]. 

Shortage costs are either of the back order type (in [$]/[Q][t]) or of 
the lost sales type (in [$]/[Q]). 

The replenishing cost is in [$] . It does not depend on the amount 
replenished . 

The reviewing cost is in [$]. It is incurred every reviewing period. 
2,k Policies 

Three inventory policies may be used. The type of policy end its 
parameters are supplied as inputs during £he execution of the program. 
Various policies and parameters may be examined for each set of demand, 
leadtime, and costs data. 

Policy 1 is the reorder point - lot size system. Policy 2 is the 
reorder point - order level system. Policy 3 is the scheduling period - 
order level system. For each policy one has to specify when inventory is 
ordered (reorder point or scheduling period) and how much is to be ordered 
(lot size or order level). 

The reviewing period is also subject to control and is the third 
parameter which is supplied for each policy. 
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2.5 Simulation 



A. The number of periods to be simulated is supplied as input during 
the execution of the program. The number of periods for which detai' ^ 'e 



to be 



B. Demands are simulated over a specified number of periods so that 
their distribution is equal to the distribution of demand supplied as data. 
The identical demands are used for each simulation. The number of periods 



is supplied as data, 



3. Presentation of Data and Inputs 



3.1 Cost Data: T9,C9,P9,H,C2,C3,C& 



A. The Period T9 



The basic time unit is the period T9« It is measured in days 
For example, T9 = 7 days. In this case, then, the period is one week. 




B. The Unit Carrying Cost : Cl 

The unit cost of the item is C9 3 e.g., C9 = $1235.71^ Per ton. 

The annual percentage of carrying cost is P9; e.g., P9 = 20$. The unit 
carrying cost. Cl, is therefore equal to T9 * C9 * P9/ 360 * 100. For example. 
Cl = 7 * 1285.71k * 20/36000 = $5 per ton per week. 



The program assumes that the year is composed of 360 days and that Cl 
is measured in [$]/ [Q] per period. 

C. The Unit Shortage Cost: C2 



Two cases are distinguished through the index H: shortages are made 
up (the back order case), H = 1; and shortages are not made up (the lost sales 
case), H=2. If H=1 then C2 is in [$]/[Q] per period, e.g. $50 per ton per week, 
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If H=2, C2 is in [$]/[Q], e.g., $50 per ton, 



D. Replenishing and Reviewing Costs: C3> C4 
The replenishing cost is C3 per replenishment; e.g., C3 = $40. 
The reviewing cost C4 is incurred every reviewing period W. E.g., if 
C4 = $2 and W = 3 weeks, then a reviewing cost of $2 will be incurred 



every 3 weeks. 



E. The Cost Data Statement 



The cost data is presented in Line 9010 : 
9010 Dm < T9,C9,P9»H,C2,C3»C4 > 



For example, 

9010 om 7, 1235.71^20, 1,50,40,0 



3*2 Leadtime Data: 1,L or 2,LL,L,L3>F(Ll) > F(L1+L3)> * * • >F(l) 

The index B distinguishes between deterministic leadtime (B=l) and 
probabilistic leadtime (B=2). When B=1 then leadtime is L periods; e.g., 
L=0 weeks. When B=2, the minimum leadtime is LI, the maximum leadtime is 



L, and the leadtime step is L3« The cumulative distribution of leadtime is 
given by the non-decreasing non-negative sequence F(Ll)» F(L1+L3)> • • • »F(L-L3)> 
F(L) where F(L) = 1. For example , if leadtime is 1,3>5> or 7 weeks with 
equal probabilities, then B=l, Ll=l, L=7> 13=2, F(l)=0.25, F(3)=0.50, F(5)=0.75> 
and F(7)=l* 



For deterministic leadtime the data statement is 

9020 Dm < 1,L > 



For example. 



9020 DA.TA. 1,0 
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For probabilistic Leadtime the data statement is 
9020 JATk < 2,U,L,I3,F(U),F(U+Io)»...,F(1) > 

For example, 

9020 DATA. 2, 1,7,2, 0.25, 0.50,0.75*1. 

3.3 Demand Data: l,X(l),X(2), . . . , or 2,X1,X2,X3>G(X1),G(X1*X3), . . . ,G(X2) 

The index C distinguishes between deterministic demand (C=l) and 
probabilistic demand (C=2) . When C=l, the demands are X(l), X(2), . . . ,X(N^). 
For example X(l)=5, X(2)=l, X(3)=k. (in this case n 4=3, but this value does 
not have to be supplied.) The demands 5»1»^>5,1,^,5, . . . are used in the 
simulation. 

When C=2, the minimum demand is SI, the maximum demand is X2, and the 
demand step is X3* The cumulative distribution of demand is 
G(X1), G(X1+X3),...,G(X2-X3), G(X2)=1. For example, if the probability 
distribution of demand is P(0)=0.05, P(2)=0.24, P(4)=0.38, P(6)=0.21, and 
P(8)=0.12, then C=2, 3CL=0, X2=8, X3=2, G(0)=.15, G(2)=.29, G(4)=.67, G(6)=.88, 
and g(8)=1. 

For deterministic demand the data statement is 
9030 DATA, < l,X(l),X(2), . . . > 

For example, 

9030 XATA, < 1,5>1>4 > 

For probabilistic demand the data statement is 
9030 DATA. < 2,X1,X2,X3,G(X1),G(X1+X3),...,G(X2) > 

For example , 

9030 DMA. 2,0,8,2, ,05, .29, .67, .88,1 
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3.4. Demands Cycle : u4 



When demand is deterministic (C=l) the demands cycle every n 4 
periods, n 4 does not have to be specified. 



When demand is probabilistic (C=2), demands also cycle ever*- n4 



periods. The program generates random demands so that their distribution 
over N4 periods will equal the given distribution G(Xl), . . . ,G(X2) . The 
number of periods w4 should be selected accordingly. It is given in the 
data statement 



9010 TATA < N4 > 



For example, 
9010 JAVA 100 



3.5 Summary of Data Statements 



The general format of the data statements is : 



9000 DMA < N4 > 

9010 Dm < T9,C9,P9>H,C2,C3,C4 > 

9020 DMA < 1,L > or 

9020 nm < 2 ,Ll,L>L 3 ,F(Ll),F(Ll+L 3 )> . • .>f(l)=i > 

9030 TATA < 1,X(1),X(2),... > or 

9030 Dm < 2,X1,X2,X3,G(X1),G(X1+X3),...,G(X2)=1 > 



For example. 



9000 Dm 100 

9010 DMA 7.1285.71^.20.1.50,1(0,0 
9020 DMA 1,0 or 

9020 nm 2,1,7,2,0.25,0.50,0.75,1 

9030 DMA 1,5, 1,4 or 

9030 DMA 2, 0,8, 2, .05, .29, .67, .88,1 




£l 













3.6 Policy Inputs: PQLICY,WIEN,HOJ MUCH, REVIEW 

During the execution of the program* the user has to supply the 
policy, J, and its parameters V(l), V(2), and V(3). The policy J can 
be 1,2, or 3. When J=1 (the reorder point-lot size system) V(l) is the 
reorder point and V(2) is the lot size. When J=2 (the reorder point -order 
level system) V(l) is again the reorder point but V(2) is the order level. When 
j=3 (the scheduling period-order level system) V(l) is the scheduling 
period and V(2) is the order level. For any policy,V(3) is the reviewing 
period: the number of basic periods in the reviewing period. 

_i:anples of inputs are : 





POLICY , WHEN , HOW 1'IUCH, REVIEW 


A. 


1 


, -2 , 


10 , 1 


B. 


2 


, 0 , 


lfc , 2 


c. 


3 


» ^ » 


10 , 



In Case A inventory is reviewed every period . When the amount on hand 
and on order is -2 or below, an order is placed for 10 units. 

In Case B inventory is reviewed every 2 weeks • When the amount on hand 

or below 

and on order is zero, /an order is placed to raise this amount to 14 units. 

In Case C the amount on hand and on order is reviewed every 4 periods. 
Orders are then placed so as to raise the amount to the level of 10 units. 



3.7 Simulation Inputs: SIM RUN > PRIMTS 

When the user supplies the policy inputs he also has to supply the 
size of the simulation run, M, and the number of periods, N, for which 
detailed information should be printed. Obviously, N ^ M. When N=0 
no detailed printouts are given, only the summary of results is printed. 

When N=l, only a brief summary is given. For example, 

SIM RUN,PRErfS 

A. 10 , 10 

B. 100 , 0 

C. 500 , -1 

In Case A, 10 periods will be simulated and details for each period 
will be printed. In Case B, 100 periods will be simulated, no details 
for individual, periods will be given, but a standard summary of results 
\rill be printed. In Case C, 500 periods will be simulated and only a 
brief summary will be printed. 

3.8 Optimization Inputs: ADV VARIABLE, STEP 

Together with the policy and simulation inputs the user has to supply 
2 optimization inputs. The first optimization input, 0, is a policy variable 
and can be )6,1,2, or 3. When Q=fi then the optimization option is ignored. 
Otherwise, three consecutive simulations are executed. The same input 
parameters are used in every simulation with the exception of the policy 
variable 0 • 
















10 



The optimization input, X3> is a step for the variable 0, In the 
three consecutive simulations the variable assumes the values 

♦ 

V(0),v(0) + X3, and V(0) + X3 + X3 respectively. This capability is 
useful for checking optimal solutions. 

Examples of optimization inputs are: 

ADV VARIABLE, STEP 

A. 0 ,0 

B. 3 .1 

C. 2 ,5 



i i 

i 




In Case A, the optimization option is ignored, and only one simulation 
will be given. (The step, too, is ignored.) In Case B, three simulations 
will be given for V(3)» V(3) + 1, and V(3) + 2. In Case C, the three 
simulations will be for V(2), V(2) + 5» and V(2) +10. 

3.9. Summary of Inputs 

The general format for request of inputs is : 



POLICY, WHEN, HOW MUCH, REVIEW, SBI RUN, PRINTS, ADV VARIABLE, STEP 



? < J ,v(l), 


V(2) 


, v(3) , M 


, N , 


0 


, X3 > 

\ \ 

\ : 

i 

V ; 

l 


Examples of inputs 


are : 






? 1 ,- 2 , 


10 


, 1 , 10 


, 10 , 


0 


• 0 (A.) j 


? 2 , 0 , 


14 


, 2 , 100 


, 0 , 


3 


. 1 (B.) j 


? 3 » 


10 


, 4 , 500 


, -1 , 


2 


• 5 (C.) : 
















r3~ar^.v??t';rr 



1 



In Case A the reorder point is -2, the lot size is 10, and the 



reviewing period is 1. There will be 10 periods in the simulation and 



the details of each of the 10 periods will be printed. No variable will 



be advanced - only one simulation will be executed. 



In Case B the reorder point is 0, the order level is 14, and the 
first reviewing period is 2. There will be 100 periods in each simula- 
tion and only the summary of results will be printed. In the second 



simulation the reviewing period will be 3, and in the third it will be 4. 



In Case C the scheduling period and the reviewing period are both 



4, and the first order level is 10. iSach simulation will have 500 periods 



and only brief summaries will be printed. In the second simulation the 



order level will be 15, and in the third it will be 20. 



. Illustrations 



Two inventory systems ape treated in this section. They illustrate 

3?Jrai 



an the capabilities of Program ISL-4 and they provide examples of actual 



computer outputs, 



i 



4.1 Illustration 1 



In an inventory system 
and have the probabilities. 



the weekly demands are in units of 2 pounds 
p(0)=.05, P(2)=.24, P(4)=.38, P(6)=.2l, and 



P(8)=.l2. Leadtime is zer 



3 . The cost of one pound is $1285 .714, and the 



ap nnai carrying cost is 2ch of this cost. Shortages can be made up, but 



I 



there is a penalty cost oj $50 per pound per week. The ordering cost is 

$1K). j 



V 








i 



Simulate the system using three different policies, 

v> 

The data, the computer output, the three inputs, and the results 
are given below: 



9000 DATA 100 
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£ 
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POLICY# WHEN# HOW MUCH# REVIEW# SIM RUM# PRINTS# ADV VARIABLE# STEP Scagple 4.13 
? 3 # 3 # 12 #2# 600 # 7 # 0 #0 



f 



l 






SCHEDULING PERIOD * 3 ORDER LEVEL » 12 REVIEWING PERIOD * 2 



PER BEG DEM END AVA. ORD REC 



CAR SHO 
ROUNDED OFF 



REP REV 



I 



I 



■> 



f 



1 

2 

3 

4 

5 

6 
7 



12 

6 

4 

12 

8 

6 

12 



6 

2 

2 

4 

2 

4 

2 



6 

4 

2 

8 

6 

2 

10 



6 

4 

2 

8 

6 

2 

10 



0 

0 

10 

0 

0 

10 

0 



0 

0 

10 

0 

0 

10 

0 



9 
5 

3 

10 
7 

4 
11 



0 

0 

0 

0 

0 

0 

0 



0 

0 

1 

0 

0 

1 

0 



0 

0 

1 

0 

0 

1 

0 



FOR A SIMULATION OF 600 PERIODS: 



CARRYING 



SHORTAGE 



REPLENISH 



g 



AVERAGES 
UNIT COSTS 
COSTS 



5*98 5 
5* 

29*925 



• 315 
50 

15*75 



•333333 

40 

13*3333 



REVIEW 

( 

•333333 



0 

0 



TOTAL COST PER PERIOD* 59.0083 



* ** * * **** ******** ************* ******* ***** ********** * ********** 



POLICY# WHEN# HOW MUCH# REVIEW# SIM RUN# PRINTS# ADV VARIABLE# STEP 
7 STOP 



RAN 



19 SEC* 



The following points should be noted: 



a. Line 9000 on Page 12 assures that 100 random demands are generated. 
These demands follow the given distribution and are used in the 
identical order in Examples 4.11, 4.12, and 4.13. 




</«/ 






15 
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mm 



b. In Line 9010 on Page 12, H=1 assures that shortages are made up 
and that the reviewing cost is zero. 

c. Since the generated demands follow precisely the given distribution 
the mean and the standard deviation on Page 12 are also the mean and 
deviation of the distribution. 

d. The abbreviations in Examples 4.11, 4.12, and 4.13 are: 

PER = period, BEG = begin, DEM = demand, END = end, AVA = available, 

ORD = order, RSC = receive, CAR = carrying, SHO = shortage, 

REP = replenishment, REV = review. 

e. In Example 4.11, the average amounts carried in the four weeks are 

3 >0,5, 2, and 0 pounds. Hence, the overall average is (3+0+5+2+0)/5 = 2 
pounds, as can be seen in the summary of the results. In a similar way 
other averages axe computed. 

f. In Example 4.12 no shortages occur during the first 8 weeks. However, 
some shortages occur during the remaining 92 weeks, so that the overall 
average is .146667 pounds. 

g. In Example 4.13 the input specifies a scheduling period of 3 weeks, and 
a le viewing period of 2 weeks. The program disregards the value of the 
reviewing period whenever the third policy (j=3) is used. It always 
assumes that in this case the reviewing period and the scheduling periods 



are one and the same. 



w «, n «-.*», - .-t. •»., - - . -■» - -• 
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4.2 Illustration 2 



In an inventory system demands cycle every 7 weeks in the following 
order: 4,2,6, 4, 8, 2, and 4 tons. Leadtime can be zero, one, or two weeks 
with respective probabilities of .2, .6, and .2. The cost of one ton 
is $1285.714, and the annual carrying cost is 20$ of this cost. Whenever 
demand cannot be satisfied, there is a loss of $50 par ton. The replenish* 
ing cost is $40 and the reviewing cost is $1.5* 



Simulate the system using three different policies. For Policy 1 
detailed printouts are desired for 15 weeks. For Policy 2 only a summary 
is desired, and for Policy 3 only a brief summary is desired. 



The data, output, input, and restats are given in Pages 17 and l8. 



The following points should be noted: 

a. The demands in all examples follow the pattern illustrated in Example 4.21, 

b. During the 10th week, in Example 4.21, the correct average amount in 
inventory is 4 2 /2x6 = 1.333. Only the rounded off value is given in 
the table. However, the correct amount is used for the summary. 

c. The main reason for the increase in shortages in Example 4.22 as compared 
with Example 4.21 is the change in the reviewing period. 

d. The abbreviations in the brief summary in Example 4.23 correspond to 



those listed earlier. TOT stands for 'total'. All the numerical values 



are costs per period. 

e. For an explanation of the 1.5 cost of reviewing in Example 4.23, see 
note g. of the previous section. 



u 
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9006 DAiA 10 

90 lu DA l A 7# 12 3 5.714# 20# 2# 63# Atit 1.6 

9020 DATA 2# U# 2# It *2» .6# 1 

9033 0 A l A \t At 2# 6# At -it 2# A 

<£/ 

rtE AD /• 



rtlW 
a'A I 1*. 



16L-4 



16:22 si 2 fttl 09/29/67 



INVENT0«</ 6/6 rEM6 61 i v J JLAT I ON# ONE PERIOD- 7 DA/6 

.J\!ir CC6l= 1 2 A 6 • 7 1 pEinCE' 11 //EArts 20 CA.-xK/InIG C06 1 /l >1 IT /PE,\l 00= 

6rlCxi A 3E6 NO l MADE JP.* 6H0 k 1A3£ COBl/J'Jll 2 50 

rcEPLEN I6-i £ NG CC6£ = 40; *£v/I£.vM3 C0oi= !• t> 

f 

LEAD i l VIE P.^OBAp lL I l / 



DEMAND A j 2 6 A i 2 A 

I>| 7 PEr\I0D6: p9£ncA3E OE'^Axl D = A. 23 571 STANDARD DEVIAlIO'l = 1.9 79 49 



* £ fs -fc 4c 4* * 5c 4c 4c # £ * £ «c jc * f: t« £ k ?: ;< .jc jc jc £ *c j-c * * k £ ^ * k, j: f? 



POLIO/# w 3E 'it NON t\.JCH # KEVlE*# 6 1 X rvJN# PkI n) 16# ADV 9Ak I ABLE# 6 1 £ P 

*> i • At lit 1# 16# Id# O# 0 




licamole ^.21 



*E vIEwING PERIOD = 1 



a EC 



C Ar\ 6H0 

kOJNQED OFF 



jnE P 



r<E \J 



bil- 



l'd 

7 

3 
id 

4 
7 

4 

6 

6 

1 

6 

1 

7 

4 

1 



v) 

0 

0 

4 

0 

u 



0 

0 

->*• 

J 



2 



j 

a 

0 

0 

2 



0 

1 

0 

1 

Id 

1 

U 

1 

u 

1 

0 

0 

1 

0 



i 






; s 



; 



!•? 



1 



H 



i i 



I 



I 



I? 












I 



R. 



I 



t 



! 



ftEr'L £ n! I OS 




• 466667 
40 

13. 6667 



k£ V I £ 

1 

1. 5 

1 • s 



TOi'AL COOi BEn» BEri 1 00- SO* 4444 

***J»****:T*1«**************V****** ?c^e fe'le |c-<« fe teAfe^e^Scie^e k4eie4e K4eie*|e4e iefele 

BOL l C > .HE >1 > S O.i S JCS* kEV l E.v* S I X kJ\l* hdilS, AOv VAk IABLfc* O'lEP 
■? 2 , 4 , 12 * 2 * 70'Jj 0 * 0 * 0 

:v£0:\0Ek POlSf = 4 Ci\G£s\ LEVEL = 12 k£v/1£a‘ I>I3 PERIOD - 2 



Oxsgiple 4.22 



f Cn A OlXJLAliCS OF 700 PErttGlKS: 



AVEk^ 3 £0 
j >i I l COM a 
COO to 



CArtrsri'l 3 

A. 3494 

5. 

21.747 



OSO.s I'A jE 



1 . 13256 

50 

59. 142 9 



i\£ PL E >1 1 OS 

• 2 5 5 5 7 1 
40 

1 1 . 5429 



n£Vt£.v 

. 5 

1.5 
. 75 



fCt'AL COOl r'Efv PEiSlOD* 93* 1327 

JC >. * H * * * * * * * k * * * * * * * * -Jc k ie fc * * * M * ■ je * * * * < jc « * •: it (s k * is «c *e * * «e k k < !« k k < '< k *: j: -is -k ’< 



r CL l Ct» /»' SE 3*3 C.i C*i J C 3 » >3 E V l £ .v > o I . v i a J >1, P,-U \i t 0> ADV VARIABLE# O'lEP Wangle 4.23 
? 3> 2> 1 2> 1 » 70> “l> 



i * 

V9 



VJ 



OCSEOjLInIj PE:\ I CO - 2 CrtOE.% LEVEL = 12 kEVIE.vl 3G PERIOD - 1 



CA;\= 22*39 5 5 03G= 41.4236 r\£P= 19. 4236 r^Ev= .75 



10T= 3 4. 5059 



Sc ^ jfc Jc i: * jc *tc jc J: jc ^c J; Jc jc :jc M jc ;jc ^c 4c Jc Jc jc *jc -Jc *jc ic ic >c jc jc < Jc -jc jc jc jc jc Jc ic -jc jc ^ jc jc Jc Jc .jc jc Jc Jc j: ic fC jc jc jc jc ^c ic ^c ie 

P0LlC/>.-iSE3»30W K-JCri* uEVIEa* Ot v i KrilMTS* ADV VARIABLE# OlEP 

? 01 OP • 



( > 



O l OB. 
rtEAOf. 




?* 



26 OEC. 
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5 . Optimization 



Comments on the construction of models and on the methods of 
fi nding optimal solutions are given in Manuals 2 and 3* In this 
section we only illustrate the use of the optimization feature of 
Program ISL-4 . 




I- 



! 

1 




1 

i 




We consider Illustration 1 of Section 4.1. Suppose one is 
interested in Policy 1 , and it is conjectured that the optimal reorder 
point is 0 and the optimal lot size is 8. The data, outputs, inputs, 
and results are as on Pages 20 and 21. 

The inputs in Example 5.1 have been selected so that the lot size 
is fixed at 8, and three reorder points are used: 0,2, and 4, the step 
being 2. One notes that the total cost for a reorder point of 2 
( 48 . 7258 ) is lower than the costs for reorder points of 0 and 4 (48.9925 

and 55.9575). Similarly, in Example 5.2 the reorder point is fixed at 2. 
The cost for a lot size of 8 ( 48 . 7258 ) is lower than the costs for lot 
sizes of 6 and 10 (49.8742 and 48.75). A further check (Example 5 . 3 ) 
indicates that for fixed reorder point of 2 and lot size of 8, a 
reviewing period of 1 week has the lowest cost ( 48 . 7258 ). It may thus 
appear that the conjecture is true. Unfortunately, this is not the 
optimal solution, as can be seen from Example 5.4. 

Incidentally, the long range expected costs can be shown to be : 



Lot Size 


6 


8 


10 


12 


Reorder Point 


-2 


81.717 


67.762 


60.100 


56.658 


0 


53.583 


48.750 


46.890 


47.317 


2 


50.017 


48.575 


48.750 


50.533 


4 


56.533 


55-962 


56.660 


58.792 



9 9 







r APh 

k£A ut . 



9000 L)A 1 A 
9013 DA 1 a 
90 P 0 dm a 
9 030 DAI A 
KEY 

KEAuY* 



100 






7# 


12 3 3* 


714# 20# 


1 # 


0 




2# 


0# 8# 


2 # .0 5# 



1 # 50# 40# 0 

.29# .67# •*«# 1 



KlM 
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i 6L-4 



16:35 v:2 FkI 09/29/67 






INv/ENfUKY SYSitlMS SIXULAIIJM# WK Pd^I iL>= 7 DAYS 
UNli C05T= 12*5*71 PEPCEM I /Y EA^r 90 CAkkYINS COST /UMl T/PICS! ! OJ)s '%• 
5H)PTA(jE$ M^L'E UP* SHORTAGE COST / ; JM I F /P£.\I0O= 50 
^PLENISHING COST= 40 TEVlEWlMG CjST= W 



X 



LtiADl l >1E= 0 



PERIODS 



OEM AMD 

0 

2 

4 

6 

8 



prCOUAUiLl TY 

.05 
• 2 4 

^ i 



21 

12 



I iM 1 O') PERIODS* AYE^AGE DEMAND = 4.22 STANDARD DEVIATION * 2*11462^ 

********* ***¥¥ ** * * ************ ********** ****** ** ****** le ** ****** 

POLICY#'-’. H EM # HO >. MUCH# PEVI E.v# SIM PUN# PP I NT S# AUV \/Ap I ABLE# STEP 



i\£0 POINT — 0 



CA* = 15.6775 SHO= 12*275 




sample $.2 



p OL I C Y * aH E<\l # HO to MUCH# REVI Ew* SIM HUM* PH I i\l TS* ADV VARI A9LE# STEP 

? 1 # 2 # 6 # 1 # 500 * - 1 * 2 # 2 



K 



I 



f 



fc 




REORDER POINT = 2 LOT SIZE = 6 REVIEWING PERIOD = 1 

CAR= 19*6553 SHO = 3*65333 rEP = 26*56 REV= 0 * TOT = 49*3742 

REORDER POI in) 1 = 2 LOT SI5E = 3 REVI E.vING PERIOD = I 

CAr= 24*7442 SHO= 2*94167 aEP= 21*04 REV= 0 TOT= 48.7258 

REORDER POINT = 2 LOT SIEE = 10 REVIEaIiNG PERIOD = 1 

CA;r= 29*67 SHO» 2*2 REP= 16*88 REv/= 0 TOT= 48.7 5 

***************** + + **-k *********** ************ + ******** *** + ♦ + *** 

POLICY* *HEN#HO.v MUCH# REVI E'.v# SIM RUN# PR I NTS* A Dv VAR I ABLE# STEP’ Zbcscrole 5.3 

? 1 # 2* 8 * 1 # 500* -1* 3* 1 * 

REORDER POINT = 2 LOT SI2E = 8 REVIEWING PERIOD = 1 

CAr= 24.7442 SHO= 2*94167 REP= 21*34 REV= 3 TOT= 48*7253 

REOruEr POINT = 2 LOT SIEE = 8 REVIEWING PERIOD = 2 

CAR= 16*7842 SHO= 28*9417 REP= 17*76 REV= 0 TOT= 63*4358 



| REORDER POINT = 2 LOT SIEE = 8 REVIEWING PERIOD = 3 

I 




p 






h 



r 



22 



******** ****:*********************** + ****************** ********* 



POLICY# /<H£iM#HO« MUCH# REMI £ri» SIN K U:M # P r\ I N T S # A U M MARI ABLE# STEP 

? 1 # 0 # 10 # 1 # 500 # ” 1 # 0 # 0 



Example 5 .4 



kEOkDEk POINT = 0 



LOT SI 5E = 10 KEMIE.vI.MG PERIOD = 1 



CAR= 20 • 41 



SHJ= 9.6 



REP = 16. R8 



REM= 0 



TOT= 46. «9 



******************************* ********:je.*4:** : ***:***.je*.* !t .*.| ;; | : * ;t ..( c ,} ; + J j : 



POLICY# WHEN# HO. 'i .MUCH# REVIEW# S I» v < RUN# PR I -MTS# ADM MARI ABLE# STEP 
? STOP 



KAN 



75 SEC. 



6. Extensions and Problems for Solution 



6.1 Extensions 



Program ISL-4 is quite a general program and allows for the analysis 
of a great variety of inventory systems. However, it is still possible 
to extend the program in many directions. Some of these extensions are 
stated in the form of exercises. 



A. The initial inventory in ISL-4 depends on the policy used. When J=l, 
the initial inventory is V(l) + V(2), otherv/ise it is V(2). Examine the 
adequacy of this initial inventory, especially when leadtime is not zero. 
Suggest other initial values to reduce the necessity for larger simulation 



runs. 



B. The occurrences of leadtime in a simulation may not necessarily follow 
the prescribed distribution of leadtime (e.g., Example 4.121). Suggest a 
more satisfactory method for generating leadtime. 

C. The cost of replenishing is independent of the amount replenished. 
Change the program so that this cost will be some function of this amount. 




/ o \s 



■ 
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]). Assume that the cost of reviewing is zero and eliminate all the 
printing statements related to reviewing. However, leave the REVIEW 
variable (V(3) in the program and allow its value to be supplied as 
input. 

E. In ISL-4 whenever an amount P is ordered, th.is amount is actually 
received after the elapse of the appropriate leadtime. Change the 
program so that the amount received is a random variable from some 
distribution one of whose parameters is P. 

F. Define a fourth policy (J=4) and allow its use during the execution 
of the program. 



6.2. Problem for Solution 

A. In a reorder point - lot size system, the amounts in inventory are 
reviewed every 2 weeks. The probability density of demand during a 
2-week period is given by f(x) » 6x(l - x), 0 £ x £ 1. The carrying 
cost is $5 per unit per week, the shortage cost is $5 per unit per week, 
and the replenishing cost is $20. Convert the density to a discrete 
distribution and show that the optimal lot size and reorder point are 
approximately 2.0 and - 0.8 respectively. 

B. Find the solution of a reorder point - lot size system when the 
reviewing period is one week, the probability distribution of demand 
during a week is P(0) = .08, P(10) = .10, P(20) = .20, P(30) = .30, 

P(40) = .16, P(50) = .10, and P(60) = .06. The carrying cost is $1 per 
unit per week, the shortage cost is $10 per unit per week, and the 
replenishing cost is $25. Show that the minimum expected total cost of 
the system is $49.35 per week. 
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C. The demand during any week is either 0 or 1 with probabilities 
0.4 and 0.6 respectively. Leadtime is 3 weeks. The carrying cost 
and the shortage cost are $1 per unit per week and $5 per unit per 
week respectively. The replenishing cost is $2. A (t,Z) policy is 
used. Find the optimal scheduling period* the optimal order level, 
and the corresponding expected minimum total cost of the system. 

D. In a reorder point - lot size system inventories are reviewed 
every week. Leadtirae is 3 weeks. During any week there is a demand 
for one unit with the probability 0.6. Otherwise there is no demand. 

The carrying cost is $1 per unit per week, the shortage cost is $5 
per unit per week, and the replenishing cost is $2. Find the optimal 
reorder point, the optimal lot size and the corresponding minimum 
expected total cost. 

E. In an inventory system with a reorder-point-order-level policy the 
reviewing period is 2 weeks. The probability distribution of demand 
during the reviewing period is P(0) = 0.20, P(0.5) • 0.24, P(l,0) = 0.40, 



p( 1 . 5 ) - o.l6. The carrying cost is $5 V er unit per week. The shortage 
cost is $10 per unit per week. The replenishing cost is $40. What is the 
optimal reorder point and what is the optimal order level? 

F. In an inventory system the reorder point is 2 and the order level is 5. 
The probability distribution of demand on any day is P(0) = 0.5, P(l) = 0.3> 



( ) 



p(2) * 0.1, and P(3) = 0.1. 

(1) Show that the probability of a replenishment on any day is 25/108. 

(2) By chan ging Program ISL-4, show that on the average there will be 
227/54 units in stock at the ^beginning of each day. 










/<w 



a * ' 






G. In an inventory system the reviewing period is 1 week. Demand 
during the reviewing period has the distribution P(0) ■ 0.2, P(l) = 0.4 
p(2) = 0.4. The carrying cost is $1 per unit per wi.ek, the shortage 
cost $1 per unit per week, and replenishing cost $1 per replenishment. 
Compare the minimum costs of the sytem for each of the three policies. 
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7. Program ISL-^ 



I SL- 4 



1110 
1 130 
1160 
1170 



OREM BY E.NADDOR 

1000DIMFC 1 1)#GC20)#HC20)# VC3) 

1 002DIMXC 500) 

1 005READN4 

1 0 1 0READT9* C9, P9, H, 02, 03, 04 
1 020L ETC 1 = T9 *C9*P9/ 36000 
1030READB 

1 040IFB*2THEN 1080 
1050REAOL 
1070GOTO 1160 
1 080READL 1*L*L3 
1100 FOR I s LI TO L STEP L3 
READ FCI) 

NEXT I 
READ C 

IF C = 1 THEN 1230 
1 180READX 1*X2#X3 
1 185READ6CX1 ) 

1 190LETHCX1 > S INT( GCX1 >*N4+»01 > 

1 195F0RI=X1+X3T0X2STEPX3 
1200READGCI) 

1205LETHCI)=INTC<G<I>-GCI-X3>)*N4+.01> 

1210 NEXT I 

1220GOTO 1350 

1230LETX 4=0 

1240LETX 5s 0 

1250LETI=1 

1260 READ XCI) 

1270 IF XCI) = 9999 THEN 1320 
1 280LETX4=X4+XC I ) 

1 29 0L ETX 5=X 5+X C I )*XC I ) 

1300 LET 1*1*1 
1310 GOTO 1260 

- 1 



1430 



1320 LET N4 = I 
1 330GOTO 1470 
1 3 SOL ETX 4s 0 
1 360LETX 5=0 
1 370FORI* 1T0N4 
1 380L£TX*X1 
1 390LETY=RNDC 0) 

1 400IFY<GCX) THEN 
1 41 0LETX=X+X3 
1 420GOTO 1400 
1 430IFHCX)*0THEN 
1 432LETXC I )=X 
1435LETHCX)=HCX)-1 
1 440LETX4=X4+X 
1450LETX5=X5+X*X 
1460NEXTI 

1 470LETX5=SQRC CX5-X4*X4/N4)/N4) 



1380 




/o i 
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ISL-4 CONTINUED 



• • t 



14S0LETX4=X4/N4 

2000PRINT 

2010PRINT”INVENTORY SYSTEMS SIMULATION# ONE PERI0D=”T9J "DAYS” 

2020PRINT 

2030PRINT M UNIT CO ST=”C95 ”PERCENT/YEAR=”P95 "CARRYING COST/UNI T/PERIOD=”Cl 

2040PRINT 

2050IFH=2THEN 2080 

2060PRINT”SHORTAGES MADE UP* SHORTAGE COST/UNI T/PERI0D=”C2 
2070GOTO 2090 

2080PRINT M SHORTAGES NOT MADE UP* SHORTAGE COST/UNI T=”C2 
2090PRINT 

2095PRINT M REPLENI SHING C0ST=”C35 "REVI EWING C0ST=”C4 

2096PRINT 

21 00IFB=2THEN 2130 

21 10PRINT M LEADTIME= M LI "PERIODS” 

2120GOTO 2190 

2 1 30PRINT”LEADTIME"# "PROBABILI TY" 

21 40PRINT 
21 50PRINTL1#F(L1 > 

21 60FORI=L 1+L3T0LSTEPL3 
2 1 70PRINTI# F< I >-F< I-L3) 

2180NEXTI 
2 1 90PRINT 
2200PRINT 
22 1 0PRINT"DEMAND 
2220IFC=2THEN 2270 
2230FORI=1TON4 
2240PRINTXC I >5 
2250NEXTI 
2260GOTO 2350 
2270PRINT”PROBABILITY' 

2280PRINT 
2290PRINTX 1# 

2300PRINTGCX1 > 

2310FORI=X1+X3TOX2STEPX3 

2320PRINTI# 

2330PRINTGC I )-GCI-X3) 

2340NEXTI 

2350PRINT 

2360PRINT”IN"N41 "PERIODS: AVERAGE DEMAND ="X 45 "STANDARD DEVIATION =”X5 
3000PRINT 

3010PRINT” 5 * :5 * ! *** 5 J ! *** 5 l' s i«******H:*4:**s|c**^***4c****^**************45*^**** 3 |e****V 

3020PRINT 

3030LETE= 1 

3040PRINT”POLICY# WHEN# HOW MUCH# REVI EVJ# SIM RUN# PRINTS# AD V VARIABLE# STEP”5 
3050PRINT 

3060INPUTJ# V< 1># V<2># VC3)#M#N#0#X3 
/4000PRINT 

40 10PRINT" - - 

4020PRINT 






I 
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I 



*■ 



I 



n 



26 



ISL-4 CONTINUED 



4580 



,/ 



4030LET2=V<1> 

4040LETQ* V< 2) 

4053LETW= VC 3> 

4500LETJ1=0 
4510LETJ2=0 
4520LETJ3=3 
4530LETJ4*0 
4540IFJ* 1THEN 
4550LETQ2=Q 
457060 TO 4610 
4580LETQ2=2+Q 
4610FORI*0TOL+1 
4620LETRCI)s0 
4630NSXTI I/ 

5000IFJ*2TKEM 504* 

5010IFJ=3THEN 50*?3 

5020PRINT"REORD?( POINT *"2l M L0T SIZE ="Q* 

5030 GO TO 5070 jj‘- 

5040PRINT M REORD/R POINT s M 2* "ORDER LEVEL = "QJ 
50 50 GO TO 5070 /' 

5060PRINT"SCHFf/JLING PERIOD = ”2J "ORDER LEVEL = ”Qi 
5070PRIN r’REVi/WIN G PERIOD s"W 
60301 FN<=0TrtFv/ 7000 
6010PRINT / 

6020PRIN T"P/» BEG DEM END AVA 
6030PRINT" REV 



ORD REC CAR SHO 



REP M I 



|M 



6040PRINT* / / 
6050PRINT/ 
7000FORK f TOM+L 
7010 LEV 01 * 02 
LV X = XCK 
If. Q2 « Q1 
f' 02 >= 0 



ROUNDED OFF” 



7020 
7030 

7100 

7 1 IF/ 
712^ -ET 
71 ^ LET 
7 4/iGO TO 
y/j0 LET 
/i60 LET 
// 170G9TO 



RC0) 

- N4*INT<<K-0.01>/N4>> 

- X 

THEN 7150 
01 > 0 THEN 7180 

11 s 0 

12 = -CQ1 ♦ Q2)/2 
7200 

II s <01 ♦ Q 2)/2 
12=0 
7200 



r 



/ 7180 LET I2«Q2t 2/<2*X) 
7190 LET 1 1 = 01 t2/C2*X) 
72001 FH= 1 THEN 7540 
72051 FQ2>=0THEN 7540 
7210LETI2=-Q2 
7220LETQ2=0 
) 7520 LET A = 0 
7530 GOTO 7550 
7540 LET A = 02 
7550 IF L = 0 THEN 7590 




/O 1 










7560 FOR I = 1 TO L 
7570 LET A = A ♦ RCI) 

7580 NEXT I 
I 7590 LET P = 0 

! 7600 LET 13 = 0 

l 7605 LET 14=0 

| 7610 IF J = 3 THEN 7680 

I 7620IFK/W=INTCK/W) THEN 7626 
| 7 62460 TO 7800 

| 7626LETI 4= ! 

I 7630 IF A>2 THEN 7800 

I 7640 IF J = 2 THEN 7690 

I 7650 LET P = P + Q 

| 7660IFA*P<=2THEN 7650 

7670GOTO 7710 

I 7680 IF K/2<>INT(K/2) THEN 7800 

I 7685 LET 14=1 

| 7690 LET P = Q - A 

7700IFP=0THEN 7800 
7710 LET 13 = 1 
I ") 7720 IF B = 1 THEN 7780 

jp 7730 LET I = LI 

f 7740 LET Y = RND(W) 

f 7750IFY<FCI)THEN 7790 

| 7760 LET I s I + L3 

! 7770 GOTO 7750 

1 7780 LET I = L 

7790 LET RCI ♦ 1) = RCI ♦ 1) + P 
7800 FOR I = 0 TO L 
7810 LET RCI) = RCI + 1) 

7820 NEXT I 

7830 LET RCL+ 1 ) - 0 

78751 FK<=LTHEN 7920 



;► 

♦ 

> 



I 



7880 LET J1 = J1 
7890 LET J2 = J2 
7900 LET J3 = J3 
7910LETJ4=J4*I4 
7920 IF K>N THEN 
7930PRINTK5 Q1»X* 
7940PRINTI 4 
7950 NEXT K 
8000LETI 1=J1/M 
801 0LETI 2= J2/M 
8020LETI 3= J3/M 
8030LETI 4= J4/M 
8040LETK 1=1 i*Cl 
8050LETK2=I2$C2 
8 060L ETK 3= 1 3* C3 
8 070L ETK 4= 1 4* C 4 
8080LETK0=K 1+K2+K3* 



11 

12 
13 



7950 

2; A5 P) RC0) J INTCI !+♦ 5)1 INTC I2+. 5)1 1 3* 



/ O 








ISL-4 CONTINUED 



8090PRINT 

8100 PRINT 

8 1 10IFN<0THEN 8260 

8 120PRINT "FOR A SIMULATION OF "MJ "PERIODS: " 
8130 PRINT 

8 1 40PRINT" "# "CARRYING"# "SHORTAGE"# "REPLENI SH"# 
8 1 50PRINT"REVI EW" 

8160PRINT 

8 1 70PRINT" AVERAGES"# 1 1 # 1 2# 1 3, 

8 180PRINTI 4 

8 1 90PRINT"UNI T COSTS”# Cl # C2# C3# 

8200PRINTC4 

82 1 0PRIN T"CO STS"# K 1 # K 2# K 3# 

8220PRINTK4 

8230PRINT 

8240PRINT"TOTAL COST PER PERIOD="K0 
8250GOTO 8290 

8260PRIN T"CAR= "K I # "SH0="K2# "REP= "K 3# 

8270PRI N T"REV="K 4# 

8280PRINT"TOT="K0 
8290IFO=0THEN 3000 
8300IFE=3THEN 3000 
831 0LETE=E+1 
8320LETV<O)=VCO)+X3 
8 330 GO TO 4000 
9000DATA100 

9010DATA7# 128 5.7 14# 20# 1# 50# 40# 0 
9020DATA1# 0 

9030DATA2# 0# 8# 2# • 05# • 29# • 67# • 88# 1 

9998DATA9999 

9999END 
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Appendix E 

FORTRAN Listings and Data 



Page 



ISL-1 : Balancing Carrying, Shortage, and 
Replenishing Costs 



Listing 

Data 



ISL-2 : The Deterministic Reorder Point- 
Lot Size System 



Listing 

Data 



ISL-3 : The Probabilistic Reorder Point- 
Lot Size System 



Listing 

Data 



112 

113 



115 

117 



119 

123 



ISL-4: A General Inventory Systems Si mu lation 

Listing 

Data 



125 

131 



111 
















I 5 L -1 ; Balancing parryin g > Shortage > and Replenishing Coots 

ISL-1 BY E*NAODOR ANO OLE BRAATEN 
FORTRAN VERSION BY RlCHARO SACHER 
THE JOHNS HOPKINS UNIVERSITY# BALTIMORE MO* 21218 

BASED ON MANUAL ISL-l* SEPTEMBER# 1967 

OIMENSION 01501*61 SOI #RI 501 #BI 501 
REAL 11* 12# 13* III* 121# 131 



INITIALIZATION 

15 OEMAND-O. 
REPIEN-O. 

1 1 - 0 . 

12 * 0 * 

13-0. 

00 20 1*1# 60 
20 R ( 1 1 * 0 * 



REAOING OF OATA IN-PERI00S#0m*06MAN0#Cl#C2#C3*UNI T COSTS! 



READ! 5# l I N 

REA0(5*2I (Dili* I •!* N I 
READ (5*21 Cl# C2# C3 



READING OF DECISIONS DATA* BMI-BEGIN* I ,R( I l-PERIOD ANO REPLENISHMENT 



REAOI 3*21 BUI 
XN ® N 

00 125 I-1#N 

125 DEMAND-DEMAND ♦0(11 
130 READ(5»Bi I# Rill 
IF (II 150*100*100 
100 REPLEN * REPLEN ♦ Rill 
GO TO 130 

150 IF (REPLEN - OEMANDI 999*200*999 

MAIN PROGRAM 

200 DO 285 I-1*N 

E ( I l»B( 1 1 ♦R ( I l-DI I I 
IF ( E( 1 1 1 260*250* 250 
240 12 • 12 - EM > 



• 12 - 
TO 270 

• II ♦ 

( R( I 1 1 

• 13 ♦ I. 



EMI 

285* 285* 280 



BUM l*E( 1 1 

AVERAGES 

Ill-Il/XN 
121* I2/XN 
I31-I3/XN 
COST l a I l l*C I 
C0ST2*I2i*C2 
C0ST3-I 31*C3 
COST ACCOST l ♦ 
COST I 1-COST 1*XN 



C0ST2 ♦ C0ST3 



112 















ISL-1 (Cont'd) 



C0ST2 l*C0ST2*XN 
COST3l*COST3«*XN 
COSTA 1 - COST 1 1 ♦ COST2I ♦ COST3I 
WRITEI6.3) 

MR|TE!6*A) (l,em,RII).0(|l,EU). !• UN) 

WRITE(6*5) 

MRITE(6*6) 1 11, I2I,I3I*CI*C2*C3 

WRITE (6* 7) COST l,COST2*COST3*COSTA*COST l 1 «COST2I *COST3I • COSTA I 
GO TO 1000 
999 WRITE(6* 10) 

C 

C NEW OAT A 

C 

1000 GO TO 15 
C 

C FORMAT STATEMENTS 

C 

1 FORMAT (81 10) 

2 FORMAT (8FI0.2) 

3 FORMAT ( 1HI,//IAX«7H PERIOD* 18X.5HBEGIN* 18X,9HSH0RTAGES, 17X, 

I 6H0EMAN0,21X, 3HEN0// ) 

A FORMAT ( IH0*6X,II2*AF25.2) 

5 FORMAT ( /////36X*9H CARRYING* I7X* 9HSH0RTAGES* 13X* IAHREPLENISHNENTS* 
I I6X*5HT0TAL// ) 

6 FORMAT ( IH0*6X* 7HAVERAGE* 5X* 3F25.2//6X* 10H UNIT COST*3X*3F25.2 ) 

7 FORMAT t IH0*6X* I IHCOST /PERIOD* IX* AF25.2/////6X* l IH TOTAL COST, 

I 2X*AF25.2) 

8 FORMAT! I10.F10.2) 

10 FORMAT ( 1HI */////// AOX, 52HREPLEN I SHMENTS ARE NOT EQUAL TO OEMANOS. 



20* 30. 20. 30. 0. 0. 

20. 20. 

10.00 





12 


















LO. 


20. 


20. 


30. 


• 

O 

CM 


30. 


0« 


0* 


M.I. 


AO. 


30. 


20. 


20. 










• 20 


5.00 


10.00 












V 


0. 


















I 


o 

• 

o 

* 















3 AO.O 

5 AO.O 

7 AO.O 

9 AO.O 

II AO.O 







12 

10. 


20. 


20. 


30. 


• 

O 

CM 


30. 


0. 


0. 


fe*. 3 


AO. 
• 20 
10.0 


30. 

5.00 


20. 

10.00 


20. 













ICHECK 


OATA. ) 


C 


END ... 




SO AT A 


12 




<4.i 


10. 

AO. 


20. 

30. 


£«. 1 


• 20 


5.00 


0. 






I 


60. 




A 


60. 




7 


60. 




10 


60. 



oQO 
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ISL-1 (Cont'd) 



1 


90.0 






3 


90.0 






5 


90.0 






7 


90.0 






9 


90.0 






II 


90.0 






t 39 _ 









12 








10. 


20. 


20. 


30. 


44,4 90. 


30. 


20. 


20. 


.20 


5.00 


10.00 




10.0 








II 


230.0 






r-93_ .... 








12 








10. 


20. 


20. 


30. 


> . 90. 


30. 


20. 


20. 


.20 


5.00 


10.00 




0. 








4 l 


30.0 






3 


50.0 






5 


50.0 






9 


70.0 






II 


90.0 






-99 








12 








$5.z l°* 


20. 


20. 


30. 


, 90. 


30. 


20. 


20. 


.20 


5.00 


10.00 




10.0 








2 


O 

• 

o 

* 






9 


80.0 






9 


70.0 






II 


50.0 






-99 








, 7 








' 5.0 


5.0 


5.0 


5.0 


/t<r ‘ 10.0 


20.0 


900.0 




0. 








1 


20.0 






5 


15.0 






-99 








7 








H* 5,0 


5.0 


5.0 


5.0 


10.0 


20.0 


900.0 





£<•7 0 * 

' I 15.0 

5 20.0 

-99 



20. 30. 



20. 30. 



20. 30. 



5.0 10.0 



5.0 10.0 
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c 

c 

c 

c 

c 

c 



IS L-2 ; The Deterministic Reorder Point--I#t Size 
I$L-2 BY E.MCCCR ANC RICHARD SACHER 
FORTRAN VERS ICN BY RICHARD SACHER 

THE JOHNS FCPK1NS UNIVERSITY, BALTIMORE, MD 21218 
BASLU CN MANUAL ISL-2, DCCLMBER , 1967 



COMMON XL, XC, XS, 11*12,13 
C I MtNS I CN CCSI5I3I, S9(3), 09(3) 
INTLGCR Xl,X2,X3,D 
REAL Jl» I 1* 12,13 









i 



\ 



A 

*> 



J 



r 



c 

C PRIMARY IRPLT 

C 

500 RE AC ( 5 , 1 ) XL ,C1 »C2,C3 

1 FORMAT ( AF 10.2 ) 

C 

C SECONDARY IRPUT 

C 

1000 REAL I 5 , A ) X $ , XC » X 1 , XU 1 ,N ,N 1 , X2 , X 3, J 1 ,0 
A FORMAT I 2F 10.2, 1 10, F 10 • 2, A 1 5, F 10.2, 1 10) 

C 

C SI MLLAT ICR ARC/CR AVERAGES AND/OR COSTS 

C 

IF (XI) 1160,1160,2000 
1160 IF 1X2) 1170,1170,3000 
1170 IF 1X3) 1100, 1180, AOOO 

C 

C TO PRIMARY CR SECONDARY INPUT 

C 

1180 IF (0 ) 5CC , 1C0C , 500 
C 

C SIMULATES TFE SYSTEM OVER A GIVEN NUMBER OF P6RI00S 

C 

2000 UR1TE(6,2) 

2 FORMAT ( IF 1 , 16X, e8HT HE DETERMINISTIC REOROER 
1PCIRT-LCT SIZE SYSTE M///) 

UR 1 Tt ( 6 , 3 ) C 1 » C2 ,C3, XU 

3 FORMAT ( 17H CARRYING COST = ,F6.2,20H PER UNIT PER PERI00,39X, 
116FSHCRTAGE COST = ,F6.2,20H PER UNIT PER PERI00/21H REPLENISHING 
2C0ST = ,F6. 2 , 18H PER REPLENISHMENT , 37X,9H0EMAN0 = *F6.2, 

3 17F UNITS PER PERIOD) 

UR I TE ( 6 , 5 ) XS , XU 

5 FORMAT ( //21X , 17H REOROER POINT = ,F6.2,32X, 1 1HL0T SIZE * *F6.2//1 
UR I IE ( 6 , 6 ) 

6 FORMAT ( IF , 50X , 19HS l M U L A T 10 N////) 

A 1 -C • 

A2 = C. 

A 3=C • 

IF (Nl) 2CAC »2C50»20A0 
20A0 UR I TC ( 6 , 16 ) 

16 FORMAT(5X,7F PER IOD , 5X , 5HBEGI N , 5X , 6H0EMAN0 , 5X , 3HEN0, 27X , 8HCARRY ING 
1, 11X,9HSFCR!AGES*6X, 13HREPLEN I SHMENT/ ) 



2050 


CC 


*090 1*1, N 




IF 


(XC1) 2059,2059,2052 


2052 


IF 


(XCl-XL) 2056,2053,2053 


2053 


11 = 


Xtil - XU/2. 




12 = 


C. 




GO 


TO 2C61 


2056 


11 


* XC 1 * *2 / ( 2 • *XU ) 



0 









I 

$ 

* 

I 

| 







115 



i 

i 



j 1 







o o o 



= ( XCl-XL )»»2/(2.»X* ) 
1C 2C6 1 



2059 



2061 



2063 



ISL-2 (Cont ’d) 

12 

CC 

I 1 = C. 

12 = Xli/2 • * XU 
XC2 * XC1 - XL 

II- (XC2 - XS) 2063,2063,2066 
I 3=1. 

Y = XC2 ♦ XC 
GC 1C 2C6S 



2066 



2069 

2070 



8 

2080 



2090 



13 = C. 

V = XC2 

IF (I - MI 207C, 2070, 2080 

14 = 13 
WRITE(6,E) 

FORMAT!/ IF 
XU = V 
A I = A I ♦ 

A2 = A2 ♦ 

A3 = A3 ♦ 



I • XCl • XU. X<)2 .11,12. 14 

, I9.F12.2, 2F10.2, 1 3X , F 20. 2 , F 19.2 , 115) 



I I 

12 



I 



COMPUTES AVERAGES FOR THE SIMULATION 



21 



FOR 



A SIMULATION OF , 14, 9H PERIODS//) 



WRITE (6,21) N 
FORMAT (/// 23H 

xn = ;j 

11 = Al/XN 

12 = A2/XN 

13 = A3/XN 

FORMAT (//67X.9F CARRY ING, 1 IX, 9HSH0RTAGES.6X , 14HREPLENI SHMENTS// ) 

FS^ATts^ C iveiiGes!llxiFil:2?2X.2F17.2//liH UNIT COSTS. 49X.FI5.2, 
1 2X.2F17.2) 

COSH = Cl* 1 1 
CCS12 = C2» 1 2 
COS 1 3 = C3» I 3 

CCS 14 = CC SI 1 ♦ C0ST2 ♦ C0ST3 
WRIlL(6iII) COSTl,COST2,COST3,COST4 
11 F0RMATI/17H CCSTS PER PER 100, 43X , F 15. 2 , 2X , 2F 1 i • 2//25H 
1PER PER ICC = ,F8 • 2 ) 

GO 10 1 16C 



10 



TOTAL COSTS 



C 

C 

C 



COMPUTES THE LCNG TERM AVERAGES OF THE SYSTEM 



17 



3000 WRITE(6,2) 

WR I T E 1 6 , 3 ) CUC2.C3.XU 
ViR I 1 1 ( 6 , 5 ) XS.XG 
WKITEIG, ID 

FORMAT! /////43X, 34H L 0 
WR I TE ( 6 ,9 ) 

CALL AVGES 
CCSU = C 1 * 1 1 
C0ST2 = C2» 12 
COST 3 = C3»I3 

CCS 14 » CCS T 1 ♦ C0ST2 ♦ C0ST3 
WRITEI6, 1C) II, 12, I3.C1.C2.C3 

CCST1,C0ST2,C0ST3,C0ST4 



ng term averages /////) 



WRITU6, 11) 
GC TO 117C 



ul6 















mu m 






1 



C 

C 

C 



CONFUTES TGIAL COSTS FOR A 
NE IGHBCHFf-CC ABCUT THtM 



ISL-2 (Cont'd) 

SPECIFIbU” KtOROER POINT AND LOT SHE ANO A 



4000 HR I I E ( 6 » 2 ) 

HR 1 1 L ( 6 , 3 ) C1,C2,C3,XU 
V> R 1 T E ( 6 , 5 ) XS,XC 
HR I IE 16. 18) 

18 FORNATI /////43X, T7H f H E 
1////56X , SF LCT SIZE//) 

CO 600 L = 1,3 
XL = L 

A9 = I XL-2. ) • J 1 
S 9 C L ) = XS ♦ A9 
600 C9IL) = XC ♦ A9 

WR1TEI6, 12) IC9IL), L « 1,3) 

12 F0RNATI1F ,23X,3F20.2) 

WR 1 I E ( 6 , 1 3 ) 

13 FGRNATI14F RECMCER POINT) 

CO 4260 1=1,3 

XS = S9( I ) 

CC 4240 J 3 1 , 3 
XC = C9|J> 

CALL AVCES 

4240 COS 1 5 1 J ) = C1«U ♦ C2«I2 ♦ C3*I3 
4260 hR 1 T E I 6 , 14) XS,C0ST5 

14 FORNATI/1F , F9 . 2 , 14X , 3F 20 . 2 ) 



TOTAL COST TABLE 



C 

C 

C 



TO FRINARY CR SECONDARY INPUT 

go tc nee 



C 

C 

c 



SUBROUTINE AVGES 

CONFUTES AVERAGES FOR THE LONG TERM AVERAGES ANO THE TOTAL COSTS OPTION 



CCNNCN XL, XC, XS, 11, 12, 13 
REAL 11,12,13 
13 = XU/XC 

IF tXS) 5C5C, 5020, 5020 
5020 11 = XS ♦ XC/2. 

12 = 0 . 

RE I L RN 

IF (XS ♦ XC) 5090,5060,5060 

11 = (XC+XS)**2/(2»*XQ) 

12 = XS««2/(2.*Xu) 

RETURN 
11 3 0 . 

12 = XC/2. “ XS 
RETURN 

JLttl 



5050 

5060 



3 



5090 



SDATA 

x« a.» . 3.0 

t.. Q 

t* *•» 5 . C 

t«. *•* 5.0 


l.C 
20 . C 
20 . C 
20 . C 


99999.0 

1 

1 

C 


36.0 

15.0 
15.0 
0.0 


1C 

100 

0 


10 

0 

0 


0 

I 

0 


0 

0 

1 


0.0 

0.0 

5.0 


0 

0 

1 


U. J.» 5.0 


99999. C 


9.0 


36.0 














f* »•« - 10 .C 


10 . C 


1 


0.0 


10 


4 


0 


0 


0.0 


0 


u *9.-10. C 


10 . C 


1 


0.0 


100 


0 


1 


1 


2.5 


1 
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ISL-2 (Coat'd) 

V- »* 5.0 i.C 9.0 0.0 

57c ”i 3.0 1C 9 0 0 0.0 0 



t».*u — | .0 


5.C 


1 


4.0 


10 


0 


l 


1 


0.5 


l 


x* »« 5.C 

t. •" -i.c 


I.C 
20. C 


9.0 

1 


_J6._0 _ 

14.0 


1C 


10 


0 


0 


0.0 


0 


- i.c 


20. C 


C 


0.0 


0 


0 


l 


l 


1.0 


l 


u 44i j.q 

-2.0 


JL*C_ 

20. C 


{ 

*0 

• 

O 

i 


36.0 

14.0 


100 


10 


l_ 


l 


0.1 


1 


Ci.id* i.c 
-2.C 


OTT 

20. C 


r.n 

i 


36.0 

4.0 


100 


8 


l 


l 


0.1 


1 


r 2400.0 
C«. AH <).C 


C • 56 
6C0.C 


99999.0 

1 


42.0 

600.0 


10 


5 


l 


l 


10. 0 


1 


r 200. C 

rr 'lM o.t 


0.04666 1 
6CC.C 


55999.0 

1 


42.0 

600.0 


10 


5 


l 


1 


10.0 


l 




_ _ 

















" 



I 

1 
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: The probabilistic Reorder Point— Ifft Siz e System 
ISl-3 HV E.MILIK AM ISUAtL PRtSSMJTO 
f- C M I RAN VERSILN BY RICHAKD SACHtR 

| HE JCFNS FCPKINS UNIVERSITY! BALTIMORE MO » 21218 



BAbtC ON PANLAL ISL-3. OCTOBER, 1967 



THIS PHCGKAP IS ECRMATTED FOR A 132 POSITION PRINTER. 

TO CHANCE FCP A 120 POSITION PRINTER, SUBTRACT 6 SPACES FROM 
E ACE FCRPAT CARC. -.*** 

CIPtNSICN rise ) ,F( 50) ,099(50) ,G( 50) ,W 160) ,V($0) .COST (3) 
INTtGER T,S,R 

REAL N(50),IC,11,I2.13,JA,L1,L2,L3,L6 
CCPPCN C,Cl,C2,C3,ll,l2,l3,0,Sl,U,N,V,X5,Wl 

3020 CC 3030 J=1,5C 
P( J )=0.C 



3030 C99(J)=C.C 
30A0 VtR 1 TE (6 , 1 ) 

1 FCHPAT( 1H1,22X,88HT HE PRO 
IPCINT-LCT SIZE S 



BABALISTIC REOROER 
Y S T E M) 



REACING ANC PRINTING OF COST DATA 



3060 RCAC ( 5f 2 ) Cl.C2.C3 

2 FONPAT ( 6F 10 .2 ) 

3070 WR 1 1 E ( 6 , 3 ) C1»C2|C3 ie< 

3 FCRPA T ( / /fcX , 16F CARRYING COST =»F9.2,20H PER UNIT PER P6RI00.35X* 

l 15FSH0RTAGE COST =,F9.2,20H PER UNIT PER PER I00//A3X , 20H REPIENIS 
2FINC CCS1 =,F9.2,18H PER REPLENISHMENT) 

3095 IRAN = 3217 
3100 Wl * 1. 



HEALING CF CEPAND OATA 

3120 RE AL ( 5 , 2 ) U.X2 
3125 JA=X2/U ♦ .Cl 
3130 J2= IFIXI J4H1 
3135 X5=C. 

3137 RE AC ( 5 . 2 ) ( P(J)» J=1,»I2 ) 
3 1 AO LO 3155 J=l,J2 
3150 XJ = J 

3155 X5=X5*U«P( J)«(XJ-1. ) 



THE EQUIVALENT DISTRIBUTION 



H70 J* J2 

1175 XJ = XJ - 1. 

1180 A = C • 

1190 B=P(J)/XJ 
1200 C99 ( J ) = A ♦ E/2. 



3210 A=A4B 
3220 J = J- 1 
3225 X J=> J- 1 . 

3230 IF (J-l) 3250,3250,3190 
3250 C99 ( 1) = PCI) ♦ A/2. 



PRINTING CF CtPANC DATA 



119 


















non non o o o o o o nonoo ooo oon 



in r " *• ^-r~^ - -* - >■'- -r - rwsr^c*. - fr<?5 ?*. <r- t W ,» vcswre ?- * ,c - x *■ nr -* kw -t?* t, 



ISL-3 (Cont'd) 



3260 hR 1 1 E I 6 , 5 ) 

5 F0RPA1 I ///A AX, 7h DEMAND, llX, 1 lHPROBABlL I TV* 10X , 1 OHCUMUL AT I VE / / ) 
3270 VI 1 ) = P( 1 ) 

3280 Gil) - CSS I 1 ) 

3290 hi 1 ) = G I 1 ) 

3300 Fill — Pill 
3310 N I 1 ) = h I 1 ) 



CCPPGTING G » h » F , N 



3320 CC 3 38C J=1,6S 

3330 FI JMPI J«ll 

3360 VIJ*1)»VIJ) ♦ FIJ+1) 

3350 G I J* 1 ) =C | J ) ♦ G99 { J ♦ 1 ) 
3360 W I J 4 1 ) = h I J ) 4 GIJ41) 

3365 IF IJ-J2) 3367,3367,3380 

3367 XJ = J 

3368 XJL = (XJ-l.l«li 

3370 WR l 1 1 1 6 » 6 ) X JL , P I J ) , F I J ) 
6 FORPATIIF ,29X,3P20.2//) 
3380 M J4 1 ) = M J) 4 W I J4 1 ) 



PAIN PRCGRAP 



6000 hRlIF. 16,98) 

98 FORPATI 1 F 1 ) 

6010 hR l IE I 6 , 7 ) 

7 FCRMT l ///6EX,35H 



//I 



HEALING CF DECISION VARIABLES, ETC (8 NUMBERS) 

IF f IRS I CF 8 NUMBERS IS LARGER 

THAN 90CC.0 TEEN NEW COST AND DEMAND DATA ARE REAO 



6060 K E A C I 5 , 8 ) S 1 , G , N V , 10 ,M , I , S,R 

8 FGRPATI2F 10 *2,110, F10«2, 6110) 

6061 IF ISI-5CC0.C) 6062 , 606 2 , 3020 

6062 X M = P 

6067 InR 1 I E I 6 , 9 ) S1,C 

9 FORP AT I///21X, 16H REORDER POINT = , F9. 2 , 62X , 10HLOT SIZE =,F9.2///1 



SIPLLAfICN ANC/CR AVERAGES AND/OR COSTS 



6050 IF I NY- 1 ) 6C6C, 5000,6060 
6060 IF IS-1) 6070,6000,6070 
6070 IF IR-1) 60CC, 7000, 6000 

SIPLLAT ICN 



5000 lnRlTE(6, 1C) 

10 FORPAT I IF ,57X,19HS I M U L A T I 0 N////) 

5080 IF IT) 5120,5120,5100 

PRINTING CF FEACINGS FOR DETAILS OF SIPULATION 



5100 UR1TEI6, 11) 

11 FORPATl 1 IX, 7F PER IOC , 5X , 5HUEG I N , 5X , 6H0EMAN0 , 5X , 3HEND , 27X , 
1 8FCARR V ING , 1 1 X , 9HSH0RT AGES , 6X , 1 3HREPLEN I SHPENT 1 

INITIALIZE S IPLLAT ICN 
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mmm 









5120 ei=IC 
5130 LI * 0. 
5140 L2 = C. 
5150 L3 = C. 
5155 L6 = C. 
5160 K = C . 



ISL-3 (Cont'd) 



V> 



FAIN SIMULA) ICN 
5170 CO 5480 J=1,M 

RANLCM CEMANC X 



5180 Y = UCRNPT ( IRAN) 

5190 J9 = 1 

5200 IF ( F(J9)-Y ) 52 10, 52 1C, 5225 

5210 J9=J9*1 

5220 GO IC 52CC 

5225 XJ9 = J9 - 1 

5230 X=XJ9»U 

5240 El=ul-X 

AVERAGES • CAR EEC ( 11) , SHORT ( 1 2 ) 

> 

5250 IF (HI) 53CC, 5300, 5260 
5260 IF (El) 533C, 5270, 5270 
5270 1 1 s ( 8 1 ♦ fc l ) / 2 • 

5280 1 2 = C . 

5290 GC TC 535C 
5300 I2=(-El-El)/2. 

5310 I 1 = C . 

5320 GO IC 5 35C 

5330 1 l = tl»»2/(2. •(El-ED) 

5340 I2 = tl»»2/(2.MB1-El) ) 

REPLENISHMENT • 13(0 OR 1) 

5350 IF (El-Sl) 5360,5360,5410 
5360 13*1. 

5370 E 2=fc 1 
5380 E2 = fc2 ♦ C 

5390 IF (E2-S1) 5380,5380,5430 
5410 I 3=C • 

5420 1 2= 1 1 

5430 IF (J-T) 5435,5435,5450 

PRINTING CF CET AILS OF SIMULATION 
5435 N I 3 = 13 

5440 WR1TE(6,12) J,B1,X, El, 11,12, NI3 

12 FORMAT ( 1FG , 6 X , 1 9, F 1 2. 2 , 2F 10 . 2 , 13X , F20. 2, F 19. 2, 1 15 1 

ACCUMULATION CF AVcRAGES 

5450 L 1 = L 1 ♦ 11 
5460 L2=L2 ♦ 12 
5470 L 3= L3 ♦ 13 
5475 L6=L6 ♦ X 
5480 Bl = E2 



SUMMARY CF SIMULATION AVERAGES 



non non 



ISL-3 (Cont'd) 



C 

5515 XM = F 
5520 1 1 = 1 1/XM 
5510 12=L2/XM 
55A0 I 3 = L 3/X 
55A5 X6=l6/XM 
55A9 WR l 1 1 ( 6 » 850 A ) M 
850A FORMAT (///5X, 20HFOR A SIMULATION OF ,I6»8H PER100S) 

5550 «RITE(6 V 85C5) 

8505 FORMAT ( / /5 A X , 7F UCMANU, ?. 3X , 8HCARRY ING , 1 1 X , 9HSH0RTAGE S, 5X , 

1 1AFREPLEMSFMENTS) 

5551 WRlTE(6t 8520 X6 , 1 1 , 12 • 13 *C 1 , C2,C 3 

8520 FORMAT! /EX, 9F A VER AGE S » 26X , 2F20 . 2 , F 19. 2 , F 1 7. 2//6X . 1 1H UNIT COSTSt 
1 AS>,F15.2»2X,2f 17.2) 

5552 C0ST1 = C1M1 

5553 CCS 1 2 = C2 M2 
555A CCSI3 = C 3 • I 3 

5555 CCS I A = CC S T 1 4 CCST2 ♦ COST 3 

5556 WR I I E ! 6 » 85AC 1 COST 1 , CO ST 2, COST 3 » COSTA 

R5A0 FORMAT! /6X, 17H COSTS PER PER 100 , A 3X ,F 15. 2 , 2X , 2F 1 7 . 2//6X , 25H TOTAL 
1C0SIS PEP PE R ICC = , F8.2 ) 

WR 11 E ( 6 ,99) 

99 FORM A T (///// 1H ) 

5560 GO 1C AC6C 

LCISC TERR AVERAGES 

6000 CALL AVER 
6030 WKlTE(6,13) 

13 FORMAT ( A8X , 3 AH LONG TERM AVERAGES) 

6035 X 6 = X 5 
60A0 HR 1 TE ( 6 , 85C5 ) 

60A 1 WR 1 1 E ( 6 , 852C ) Xfc, 1 1 , 12, 13,C1»C2»C3 
60A2 C0ST1 = C1M1 
60A 3 CCS 1 2 = C2 • 12 
CCST3=C3» I 3 

60A A COS I A = CC S T 1 4 COST 2 + C0ST3 
60A5 WR 1 1 E ( 6 , 65AC ) COST 1 , COS T 2 , COST 3 , COSTA 
WRI IE ( 6 ,99 ) 

6070 GO IC AC7C 

EXPECTEC T01AL COST TABLE 

7000 WR1TE16, 1A) 

1A FORMAT ( A3X, A7H EXPECTEO TOTAL COST TABL E//// 
1 62 X , 9H LCT SIZE) 

7035 C11=C-L 

7036 C12=C 

7037 C13=C ♦ L 

70A0 WR|TL(6,15) C11»U12,G13 

15 FORMAT! 1F0,29X, 3F20.2/21X, 1AH REOROER POINT) 

7060 S 1 = S 1-2 • »L 
7070 C = C ♦ L 
7080 CO 7180 1=1,3 

7090 S1=S1 ♦ L 
7110 C=C-3.*L 
7120 CO 7150 J= 1 , 3 
7130 C=C ♦ U 
7 1 AO CALL AVER 



noon o 



7150 COS I ( J ) = C ISL-3 (Cont 

7180 WRIIE(6, 16) SltCOSr 
16 FCRMTUFC, 9X,4F20.2) 

7390 GO 1C 4CCC 
C 

ENt 

SUBROUTINE aver 
C lMt-NSICN V ( 50 ) 

RLAL N( 50, 1C, I 1, 12, I 3* J4 , L 1 , L 2 , L 3 , L 6 
INItGER 1,S,R 

COMNCN C,C1,C2,C3,U,12,13,U,S1,U,N,V,X5,W1 
8000 B — (SI ♦ C - U/U + .01 
8005 NB = INT(B) ♦ 1 
8010 IF (NO - 1) e02C, 8050,8050 
8020 I 1-C . 

8030 GO (C 8150 

8050 IF (SI) 8100,8100,0060 

8060 A = (SI - U)/U ♦ .01 

8069 NA = 1NT(A) ♦ 1 

8080 ll=( (U**2)/C)*(N(NB)-N(NA)) 

8090 GO TO 8150 

8100 1 1= ( ( U • • 2 )/C ) »N ( MB ) 

8150 12=11 ♦ X5/2. - ( (Q+Ul/2. ) - SI 

8151 IF (12 - .0CC001) 8153,8155,8155 
8153 12 = 0. 

8155 C =(C-U)/L 4 .Cl 
8160 NO = INT(C) ♦ 1 
8200 13 = ( 1.- (L/C)»V(NU))/Wl 
8210 C-CI«ll ♦ C2 • 12 ♦ C3» I 3 
RETURN 

ENC 

FUNCTION LCRNRT(J) 

THE FUNCTION UCRNR T IS A UNIFORMLY DISTRIBUTED RANOOM NUMBER GENERATOR 
WRITTEN BY MANCELL BLLLMORE, THE JOHNS HOPKINS UNIVERSITY 
THIS FUNCTION WAS WRITTEN SPECIFICALLY FOR THE IBM 7094. THE ARGUMENT, 
IRAN, POST BE INITIALIZED TO AN 000 INTEGER NOT DIVISIBLE BY 5 AND MUST 
LIE BETWEEN C AND 4194303 



JRAN = J 

JRAN = JR AN • 205 1 

JRAN = JRAN - ( JRAN/4194304)*4194304 

J = JRAN 

U = JRAN 

Y = U/41943C3.C 

UORNRT = Y 

RETURN 

END 

SDATA 



5. 


5C . 


40. 










2. 


e. 












. .05 


.24 


• 

Ul 

CD 


.21 


.12 






-4.C 


14. C 


1 


id. o 


10 


10 


0 


C,.4n- -4.C 


14. C 


1 


10.0 


500 


10 


1 


2.0 


8.C 


c 


0.0 


0 


0 


0 


0>c 


10. C 


1 


4.0 


1000 


0 


1 


9999. C 


O.C 


c 


0.0 


0 


0 


0 


l.C 


9.C 


36.0 










^ 5.C 


5.0 












. o.c 


1 l.C 












o.c 


20. C 


0 


0.0 


0 


0 


0 


9999. C 


O.C 


C 


0.0 


0 


0 


0 
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ISL- 4 : A General Inventory Systems Simulation 
c isl-a by t .rr n r n n * ^ 

C FCRIRAN VERSION BY KICHAKO KE SSI NGER ANO JON WESTON 

C THE JOHNS HCPKINS UNIVERSITY, BALTIMORE MO. 21210 

C 

C BASIC CN MANUAL ISL-A, CCTOBEK, 1967 

C 

IN I EGER NA, 19, IH,B*L*Ll*L3*H,Xl,X2»X3,C*X*E»W,JfV»M*N»0 
INTEGER X6,X7,IA,15,1X 

REAL C9 , P9 , C2 *C 3 ,CA , C l * FNA ,l»* F * XA » XB ,XA , X5* XD *P, A* 01 *02 ,Q* R*l 
REAL 11, 12, 13, IA,J1,J2, J 3 * JA , K 1 , K2 , K3 , KA , KO 
C 

c SLIBKOUIINC HHN I I ) GENERATES RANDOM NUMBERS 

C IN THE IMLPVAL 0 TO 1. SUPPLIED BY THE USER 

C 

LETY = RRN I 12 ) 



C 

1000 

1002 

C 

c 

c 

1005 

1007 

1 

1010 

2 

1015 

1020 

1030 

3 

10A0 

1050 

1070 

1080 

A 



1110 

1120 

5 

1130 

1160 

1170 

1180 



1185 

1187 

1190 

1193 

1195 

1200 

1203 

1205 

1210 

1220 

1230 

12A0 

1250 



DIMENSION F(11)*G(20)*V(3),IH( 10) ,R(20) , IRI20) 
CIMENSICN X ( 50C ) 

RtALING CF LATA AND INITIALIZATION 



REAL (5,1) NA 

IF (NA.EC .9959 )GC TO 9998 

FORMAT (110) 

REAL (5,2)T5,C9,P9,H,C2,C3,CA 

FORMAT ( I 10, F 10. A, F 10. A, I 10, F 10 . A , F 1 0 . A , F 10. AJ 
F T9 = T9 

Cl * (FTS»C5»P5)/36000. 

REAL (5,3) e 
FORMAT (110) 

IF (B.EC.2) GC TO 1080 

REAC (5,3) L 

GO TO 116C 

REAC (5, A) L 1 , L , L 3 

FORMAT (3110 

L A = L 1 ♦ 1 

L5=L*1 

CO 1130 I=LA,L5,L3 
REAC ( 5,5 > F ( 1 ) 

FORMAT (F10.A) 

COM INGE 
R E AC I 5 , 3 ) C 

IF (C.EC.l) GC TO 1230 
REAL (5, A) X1,X2,X3 
X6*X !♦ 1 
X7=X2* 1 

REAC (5,5) C ( Xfc ) 

FNA = NA 

I H ( X6 ) = G ( Xt ) »FNA 4.01 
I A = X6 ♦ X 3 
CC 1210 I s I A , X7, X3 
REAC (5,5 ) G( 1 ) 

It) = I — X 3 

1 H ( 1 ) = (Gin — G I IB ) ) »FNA 4.01 
CCM INUE 
GC TO 1 35C 
X A = C 



X5 = C 

1=1 



O 

=ERLC 



•j 
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m. 






m. 



U60 «tn< 5 .t>Mi. gki tea> a i 

6 FORMAT (I1C) 

1270 IF ( X C 1) .EC.9999) GO fO 1320 
1275 XU = X (1) 

1280 X A = XA OC 
1290 X5 - X5 ♦ XC«XU 
1300 1= I ♦ l 

1310 GO 1C 126C 
1320 NA = I-I 
1325 FNA -NA 
1330 GO 10 1 A7C 
C 

C GENERATION CF NA RANUUM NUMBERS 

C 

1350 XA =0 
1360 X5=C 

1370 CO 1 A60 1=1. ISA 
1380 IX-XU1 
1390 Y = «RMC) 

1 AOO IF (Y.LT.G(IX)) GO TO 1A30 
1 A 10 IX = IX ♦ X 2 
IA20 GO 1C 1 ACO 

1 A 30 IF llh(lX).EC.C) GO TO 1380 
l A32 XUIMX-1 
l '*3 A HX = X ( I ) 
l A 35 I H ( IX) -IH IX) -1 
1 A AO XA = XA*FX 
1 A50 X5 = X5 ♦ FX«FX 
1A60 C0MINC6 
C 

C MEAN ANC STANDARD dlviation 

c 

1 A 70 X5 = SCRT ( ( X5-XA»XA/FNA )/FNA ) 

IA 80 XA = XA/FNA 
C 

C PRINTING CF FEACINGS AND DATA 

C 

2000 WRITE (6.7) 

7 FORMAT (/) 

2010 IriR I I E (6.8) T9 

8 FORMAT ( A2H 1 INVENTORY SYSTEM SIMULATION. ONE PERIOO *,I10,5H OAYS) 
2020 WRITE (6,7) 

2030 WRITE (6,9) C9.P9.C1 

9 FORMAT ( l 3H LNIT COST = .F10.A.16H PERCENT/YEAR = .F10-A.29H CARRYI 
CNG COST/CNIT/PER IOD = .F10.A) 

20A0 WRITE (6,7) 

2050 IF (H.EC.2) GO TO 2000 
2060 WRITE (6,10) C2 

10 FORMAT ( 50F SHORTAGES MADE UP. SHORTAGE COST/UN I T /PER 100 * ,FiO 
C • A ) 

2070 GO TO 2C5C 
2080 WRITE (6,11) C2 

11 FORMAT ( 52F SHORTAGES NOT MADE UP. SHORTAGE COST/UN I T/PER 100* , F 10 
C • A ) 

2090 WRITE (6,7) 

2095 WRITE (6,12) C3.CA 

12 FORMAT (22H REPLENISHING COSTS = ,F10.A,19H REVIEWING COSTS * ,F10 
C.A) 

2096 WRITE (6,7) 

2100 IF (B.EC.2) GO TO 2130 
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i 



t 




I 



i 

i 



i 



2110 

13 

2120 
21 30 

14 
2140 
2150 

15 



HR I It* (6,13) L 
FORPAT (12H LE AIL T IMt 
GO 1C 219C 
wR I IL (6,14) 

FGRPAT ( 2 IF 
HHIlt ( fc , 7 ) 

UK ! It (6, 15) 11, F ( L 1 ♦ 1 ) 
FORPAI (IF , I 1C,5X,F 10.4) 



isl-4 

= , 1 10 , 8H 



(Ccmt'd) 

PERIODS) 



LtAOMMF PROBABILITY) 



2155 IC = L l ♦ L 3 ♦ 1 
2160 CC 2 1 RC 1=IC,L5,L) 

2165 10 = I-L2 

2167 XA =F ( I )-F ( II ) 

2168 I AA = I- 1 

2170 v»R I IE (6,15) IAA, XA 
2180 COMINLC 
2190 HRllE (6,7) 

2200 HR I It (6,7) 

2205 IF (C.tC.2) GO TO 2270 
2210 HR 1 1 6 (6,16) 

16 FORPAT ( 7F LEPAND) 

2215 HKII6 (6,7) 

2230 CC 2250 1=1, N4 
2240 HR lit (6,17) X ( I ) 

17 FORPAT (IF ,7111) 

2250 CCMlNUt 

2260 GO 1C 235C 
22 70 HR I Tt ( 6 , 18 ) 

18 FORPAT ( 29F DEPANO PROBABILITY) 

2280 HRITl (6,7) 

2290 HR 1 1 6 (6,19) Xl,G(X6) 

19 FORPAT (IF , I IC,5X,F 10.4) 

2300 IG = X6 ♦ X 3 

2310 CC 2340 I=IG,X7,X3 
2320 1 J= I- X 3 

2325 XR = G ( I 1— G ( IJ) 

2326 I AB* I— 1 
2330 HR IIE(6, 19) IAb,Xfl 
2340 CON I INGE 
2350 HR I 1 E (6,7) 

2360 HRITE (6,20) N4,X4,X5 

20 FORPAT ( 4F IN ,110, 27H PERIODS- AVERAGE DEPANO * 
CRD DEVIATION = ,F10.4) 



, F 10.4 , 22H STANOA 



C 

C PAIN PRCGRAP 

C 

3000 HRITE (6,7) 

3010 HR 1 1 6 (6,21) m ..... 

21 FORPAT (76F ******************************** 



3020 HRITt (6,7) 

3030 6=1 

C READING CF CONTROLLABLE VARIABLES, ETC. (8 NUMBERS) 

C 

3050 r<6AC(5,23) J , V ( 1 ) , V ( 2 ) , V ( 3 ) , M,N, 0, X3 
23 FORPAT miC) 

3051 IF ( J. EC. 5995 )GC TO 1005 

3 °22 FORPAT (76F1 POLICY, WHEN, HOW MUCH, REVIEW, SIM RUN, PRINT 
CS, ADV VARIABLE, STEP ) 



12V 





I 

I 

I 

■j 

if 



1 

| 

I 

| 

I 

j 

% 

n 

S 

| 

I 

■: 










mm 



'» ,i~ t 



1 



3045 

3055 

24 
C 

C 

C 

4000 

4010 

30 

4020 

4030 

4040 

4050 

4500 

4510 

4520 

4530 

4540 

4550 

4570 

4580 

4610 

4615 

4620 

4630 

C 

C 

C 

5000 

5010 

5020 

25 
5030 
5040 

26 
5050 
5060 

27 
5070 
5072 

28 
C 

C 

C 

6000 

6010 

6020 

31 

6030 

32 

6050 

C 

c 

c 

7000 

7001 
G 



ISL-4 (Coat'd) 

WRITE 16.71 i L 

WRITE <6.24) J.VI 1 I . V ( 2 ),V( 3),M, N.U.X3 
FORMAT IIP ,819) 



INITIAL 12 AT lUN OF SIMULATION 



WRITE 16.7) 

WRITE (6,30) 

FORMAT (76H ------------- 

C ----- ) 

WRITE 16.7) 

Z- V f 1) 

C* VI 2) 

W*VI 3) 

J1=0. 

J2=0 • 

J3=0. 

J4 = 0 • 

IF IJ.EU.l) GO TO 4580 
C2=C 

GC TO 4610 

C2=Z4Q 

IK=L*2 

CC 4630 1*1, IK 
R I I ) =04. 

CONTINUE 

PRINTING OF POLICY AND ITS PARAMETERS 

IF IJ.EC.2) GO TO 5040 
IF IJ.EQ.3) GO TU 5060 
WRITE (6,25) 7,0 

FORMAT I 17H REORDER POINT = .F1C.0.12H LOT SIZE = .F10.0) 

GC TO 5070 
WRITE 16,26)2,0 

FORMAT I 17H REORDER POINT = ,F10.0,15H OROER LEVEL - ,F10.0> 

GC TO 5070 
WRITE I 6,27 )Z,Q 

FORMAT I21H SCHEDULING PERIOD = ,F10.0,15H OROER LEVEL = ,F10.0) 
WRITE 16,7) 

WRITE(6,28) W 

FORMAT I20H REVIEWING PERIOD = ,110) 



DET AILED SIMULATION HEADING 



IF IN.LE.O) GO 
WRITE 16,7) 
WRITE 16,31) 
FORMAT I68H 


TO 7COO 
PER BEG 


OEM 


ENO 


AVA 


ORD 


REC 


CAR SH 


CO REP REV) 

WRITE (6,32) 
FORMAT I59H 
C CFF) 

WRITE (6,7) 














ROUNDED 



MAIN SIMULATION 



MPLUSL*M4L 
DC 7950 K* 1 .MPLUSL 



[ 

i 

f 

I 

I 



b 

* 

I 









r 



tWBpWwiWwaitga^ 















c 

c 

70m 

c 

c 

c 

701? 

70l<* 

7015 

7016 
7020 
7030 

C 

c 

c 

7100 

7110 

7120 

71210 

7140 

7150 

7160 

7170 

7180 

7190 

72C0 

7205 

7210 

7220 

C 

c 

c 

7520 

7530 

7540 

7550 

7560 

7570 

7575 

C 

C 

c 

c 

7590 

7600 

7605 

7610 

7620 

7624 

7626 

7630 

7640 

7650 

7660 

7670 

7680 

7681 
7685 
7690 

7692 



IGL-4 (Cont'd) 

Cl AND C2* BEGIN ANC F.NL imNTCKltS 

Ci = C2*« ( l 1 

CEP<-'N0= EX 

FK = K 
FIS4 = N4 

INTFK=IMUFK-0.01)/FN41 

I XSUB-K-N4* INTFK 
FX = X I IXSUB) 

C2 = C l-FX 

AVERAGE CARRIED 3 1 It AVERAGE SHORT = 12 

IF IC2.GE.0.0) GO TC 7150 
IF ICi.GT.0.0 ) GC TC 7180 

II = 0. 

1 2-- (01+C21/2.0 
GC TC 7200 
I 1® ( CI+C2 I / 2«0 
12=0. 

GC TC 7 2C0 

l?=C2**2/ ( 2 • *FX 1 

I l=C 1**2/ I 2 «*FX 1 

1FIH.EC.1) GC TG 7540 

IF IC2.GE.0. ) GC TC 7540 

1 2--C2 

C2=0.0 

A= AFCUNT CN HAND ANC ON ORDER 
A=0. 

GC TC 7550 
A=C2 

IF ( L.EG .0 1 GC TC 7590 
CC 7575 1*2. L5 
A-A«R( I ) 

CCNT INUE 

I3 = REPLEMSM0 CK 1)14=KEVIEW (C OR l) 
Re APCUNT REPLENISHEC 

P = C. 

13=0. 

14=0. 

IF I J.EQ.31 GO TO 7680 
IFIK-W*IK/W).EQ.O) GO TO 7626 
GC TO 7800 
14=1. 

IFIA.GT.Z) GO TO 7e00 
IFU.EC.2) GC TC 7690 
P=P*G 

I F I A+P. LE .Z ) GC TC 7650 
GC TC 7710 
!Z= INT I £♦• 5 ) 

IF|K-IZ«(K/IZ).NE.C) GO TO 7800 
14=1. 

P = G- A 
IP=P 
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) 



ISL-4 (Cont'd) 

7700 IF I IP-EC -0 ) GO TO TPOrt 
7710 13=1-0 

7720 IMB.E0.1) GO TO 7780 
7730 1*L1->1 

c 

C R ANCCN LEADTIME 

C 

77A0 Y«RRN I 0) 

7750 IFlY-LT.mn GO TC 7790 
7760 I* l*L 3 
7770 GC TO 7750 

C UPCATE AMOUNTS ON CRDER Kill 

C 

7780 1 - L ♦ 1 

7790 R 1 U1) = RJ l + D + P 
7800 CC 7825 1*1. L5 
7820 Rl 1 1*R1 I ♦ i 1 
7825 CONTINUE 
7830 RIL*2)=0. 

7875 IFIK.LE.L) GO TO 7520 

C 

C ACCUMULATE 11,12,13 

C 

7080 Jl=Jl*ll 
7890 J2=J2*I2 
7900 J3=J3*I3 
7910 JA=JA*IA 
7920 IF1K.GT.N) GOTO 795C 



C 

c DETAILS OF SIMULATION 

C 

7922 IX=INT| FX4.5) 

7923 IR11 ) = 1NT|RI 1I+.5) 

792A * P* 1NT 1 P + -5 I 

7925 H 1= INT ( 1 1 ♦ • 5 J 

7926 1 1 2* INT 1 1 • 5 > 

7927 ICl=KEShESI01» 

7928 1C2=KESWES ( C2 ) 

7929 IA=KESWES<A) 

7930 I 1 3=KESWES ( 131 

7931 IIA=KESWESIIA> 

7932 KR1TEI6,79A0 )K, 101 , IX, I02r I A, IP, 
79A0 FORMAT (1116) 

7950 CONTINUE 

7999 FM=M 

8000 1 1 = J 1/FM 
8010 I2=J2/FM 
8020 13=J3/FM 
8030 IA* JA/FM 
80A0 Kl= 1 1 *C 1 
8050 K2= 1 2*C 2 
8060 K3» I 3*C 3 
8070 KA=IA*CA 

8080 K0*K1 +K24K34KA 
6090 WR1TE(6,7) 

8100 WRITEI6.7) 

0110 IFIN.LT. 01 GO TO 3260 



IR( 1), 111,112, II3,IIA 



C 

C SIMULATION SUMMARY 
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( ) 



I 



“5F 

ERIC 



8 140 
8145 
8160 
8170 
0175 
8190 
8195 
8210 
8215 
8230 
8240 
8245 
8250 



shortage 



REPLENISH 



IGL-4 (Coat’d) 

8 l ,!° Mi»‘.ti!!rc^c5 A S l MUL AT 1 ON OF . ! 10, 8HPER .005 . 

8130 wKlTh(6.7) 

WRIltiC, 61451 
FGHNATI 14X»5CPCARRY ING 
k R 1 1 E I 6 t 7 ) 

MM I f fc C 6.6175 I 1 1» I2» M > IA 

FORMAT! IX , 1 3FAVERAGES 
V,RI 16(6,6195 ) C1»C2*C3#C4 
FOHN AT I IX » 13FLM I COSTS 
WRIIt(6,8215) Kl,K2,K3»K4 
F C R i v A T I IX , 13FCCS TS 
WR I TE I 6 » 7 ) 

FORMltlXf22FTCTAL COST PER PER 1 00 = • F 1 2 • 4 ) 

GC TC 829C 



REVIEW) 



,F12.4tF12.4»F12.4»Fl2.4) 

,Fl2.4,Fi2.4 t Fl2.4,F12.4) 

f F12.4,F12.4,Fl2.4,Fl2.4) 



C 

C 

C 



BRIEF SIMAATICIS SUMMARY 



8260 

8265 



roi;Mu5!tF^^:FK!A:AHS»S=.FlO.A.AH R £P..FlO.A.AHKeV-.FI0.4, 

C4HICT=»F1C.4) 



ADVANCING a VARIABLE 



8290 IFIC.EC.O) 
8300 IFIE.6C.3) 
8310 E=£*l 
8320 VIC ) = V(C ) + X3 
8330 GO TC 4CCC 
9998 S TC P 



GC TC 
GC T C 



3000 

3000 



9999 tNC 

FUNtTICN KESWESIX) 
IF(X.GE.C.C) GC TO 
40 KEShES= INTI X-.5 ) 
RETCRN 

30 KESNES=INTIX4.5 ) 
RETLRN 

r (El P 


30 












hlNL 

IIBLCR RRN* 
STEXT RRN* 












KKN * UU vju 
RRN* 000l 
RRN«0002 


SCC 1 C T RRN • 
tiWFNC RRN • 












RRN*0004 

RRN»0005 


^ U L 1' w * * * 

(DATA 

ICC 

HI. 4.1 \ 1285 - ? “ 


20. 


1 


50. 


40. 


0.0 




C 

2 

C 8 

• 05C 

• 29C 

.67 

.88 

% 


2 












i - 2 

Ll<. 2 0 


8 

1C 


1 

1 


5 

100 


5 

8 


0 

0 


C 

c 
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